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On the basis of the Blaizot-Iancu equations, which are a local formulation of the hard thermal loop 
(HTL) equations of motion for soft fluctuating quark and gluon fields and their induced sources, 
the coupled kinetic equations for plasminos and plasmons are obtained. The equality of matrix 
elements for nonlinear scattering of a plasmino by hard particles in covariant and temporal gauges 
is established by using effective Ward identities. The model problem of the interaction of two 
infinitely narrow packets with fermion and boson quantum numbers is considered. The kinematical 
relations between wave vectors of the plasmino and plasmon are derived, when the effective pumping 
over of the plasma excitation energy from the fermion branch of plasma excitations to the boson 
branch and vice versa occur. The expression for the nonlinear Landau damping rate of a plasmino at 
rest is found, and a comparison with a plasmino damping constant obtained within the framework 
of the hard thermal loop approximation is made. The nonlinear Landau damping rate for normal 
quark excitations is shown to diverge like l/\/q^ near the light cone where q is a four-momentum 
of excitations, and the improved Blaizot-Iancu equations removing this divergence are proposed. 
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I. INTRODUCTION 



This work deals with the study of the interactions between soft-fermion and gluon excitations in an off-equilibrium 
quark-gluon plasma (QGP) within the hard thermal loop effective theory. The specific problem of the damping of 
soft, collective, fermion excitations is ultimately considered. This problem has already been addressed in the literature 
(see the review below) in the formalism of thermal perturbation theory, but nevertheless it would be interesting to see 
how it emerges in a real-time formalism based on kinetic equations for soft excitations. In addition, the calculation 
of the damping rate for moving soft excitations is still incomplete. 

The theoretical investigation of properties of plasma excitations bearing fermion quantum numbers in a quark- 
gluon plasma originated in the pioneering works by Klimov and Weldon M. The gauge-independent dispersion law 
of fermion excitations was first derived in a high-temperature region within the framework of the imaginary time 
formalism (Matsubara technique) by Klimov and within the framework of the real-time formalism (Keldysh- Swinger 
technique) by Weldon at one-loop order. The quark propagator in the case of nonzero temperature and massless 
fermions was shown to have two poles corresponding to two different dispersion relations, both with a positive 
energy. The first branch describes normal-particle excitations with the relation between chirality and hclicity at zero 
temperature. The second branch is collective excitations, where the usual relation between chirality and helicity is 
flipped 0. It has been called a plasmino Q to emphasize that, like the plasmon mode of the gluons, it is a purely 
collective branch of plasma excitations that have no analog at zero temperature. 

The damping rate of plasma modes is one of the most important characteristics of high-temperature plasma. The 
calculation of the damping rate of fermion excitations in the rigid one-loop approximation was made by Kapusta 
and Toimela They showed that the standard one- loop calculations were incomplete and resulted in an explicit 
dependence of the damping rate on the gauge-fixing condition and were even negative in some gauges. Overcoming 
this problem both in the case of fermion and boson modes is known to lead to the development of an effective 
perturbation theory put forth by Braaten and Pisarski || and by Frenkel and Taylor j(| that enables one to resum 
in a systematic way high-loop diagrams contributing to damping rates at the leading order in the coupling constant 
g and gives formal proof that resummation produces gauge-invariant results for the damping rates of both quarks 
and gluons. The calculation of the damping rate of a heavy fermion H by Pisarski is one of the first examples of the 
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application of effective perturbation theory. 

The damping rates depend in a nontrivial way on how fast the quasiparticle is moving through the thermal medium. 
Within the framework of the imaginary-time formalism using a dispersion relation method, Braaten and Pisarski || 
and independently Kobes et al. |9[| present by somewhat distinct approaches a complete calculation of the quark 
damping rate at zero momentum at leading order in g by evaluating one-loop diagrams constructed out of the 
effective propagators and vertices. It is shown that the damping rate of quark modes at rest is some pure number 
times g 2 T ', where T is a temperature. 

Much interest has been taken in the so-called damping rate of moving (hard or soft) quarks in hot QCD plasma 1 . 
The hard thermal loop (HTL) resummation renders the infrared finite damping rate of plasma excitations at rest. 
However, the damping rate of fast moving collective excitations is infrared divergent, arising from the exchange 
of quasistatic, magnetic gluons, reflecting the breakdown of the Braaten-Pisarski resummation scheme in deriving 
the physical quantities that are sensitive to 0(g 2 T) corrections. In hot QCD plasma this problem is avoided by 
the introduction of the magnetic screening mass of order g 2 T ', which is expected to come from soft-gluon mutual 
interactions pi] . An accurate calculation of the damping rate for soft moving quarks and gluons has been done by 
Pisarski in Rcf. Jl^] . He found that the infrared-singular contributions were proportional to the group velocity of the 
respective quark and gluon modes times g 2 T\a{\/g). In Ref. 1 1 3 Flechsig et al. presented a unified treatment of all 
next-to-leading-order corrections to fermion and gluon dispersion laws in QGP's, which are infrared singular due to 
mass-shell singularities. They reproduced the result of Pisarski jl2| for moving quasiparticle excitations. 

The need for further resummation beyond the scheme of Braaten and Pisarski also arises in an attempt to calculate 
the production rate of soft real photons from equilibrium QGP's iQ. This rate was shown to be logarithmically diver- 
gent owing to the so-called collinear (mass) singularities of a HTL diagram. The reason for the collinear singularities 
is known to be the use of bare (massless) hard-quark propagators that are on the mass shell. In a somewhat different 
context a similar difficulty was encountered in the research of the behavior of the next-to-leading order longitudinal 
component of the polarization function of the hot gluon system near the light cone |l5] . Here, the singularity is even 
more strong and generated by a premature restriction to soft loop momenta (the detailed study of the light cone 
singularity in the simpler gauge theory of scalar electrodynamics was made by Kraemmer et al. [ |16| ) . 

The way to solve the problems connected with collinear singularities in hot QCD plasma was proposed by Flechsig 
and Rebhan Jl7|] . They showed that the inclusion of asymptotic thermal masses for hard transverse gluons 2 and hard 
quarks removed the collinear singularities of the HTL amplitudes without spoiling gauge invariance. The examples of 
employment of the effective hard-quark propagator and effective photon-(hard)-quark vertex for deriving the finite 
soft-photon production rate from the equilibrium QGP can be found in Refs. |l7| pj| . 

Here we would like to represent a somewhat different view about the study of the above-mentioned problems, based 
solely on the kinetic theory. More comprehensively, our purpose is to construct kinetic equations that describe the 
evolution of off-equilibrium deviations in the number densities of the soft excitations (fermions and gluons), and to 
obtain damping rates from the collision terms. As we have shown in the case of soft-boson modes ^O], pi] , there is 
no direct connection between the results obtained within the framework of kinetic theory of soft excitations and the 
results derived within the framework of effective perturbation theory. Here, we would like to extend the analysis, 
carried out in Refs. on the fermion degree of freedom of plasma excitations in hot QCD plasma. Our approach 

is based on the fundamental system of equations derived by Blaizot and Iancu |2^| . These equations are obtained 
on the basis of a truncation of the Schwinger-Dyson hierarchy. They isolate consistently the dominant terms in 
the coupling constant in a set of equations that describe the response of the plasma to weak and slowly varying 
disturbances, and encompass all HTL's. This system of equations is determined on the space-time scale (gT)^ 1 . We 
show how to obtain the information connected with the damping of soft excitation QGP's by using this system, i.e., 
with the processes determined on the other space-time scale. 

In order to allow for the collective excitations with arbitrary quantum numbers, Blaizot and Iancu considered both 
fermion and boson soft fields. The pure gauge sector (when the influence of the soft-fermion field is neglected) of 
this system of equations has been studied in detail. The sector connected with a soft-fermion degree of freedom of 
plasma excitations is not actually studied. In the present paper we have shown that this sector of plasma excitations 
contains some extended information on the dynamics not only of soft-fermion modes, which are especially interesting 
for us, but also soft-boson modes interacting with them. Deriving a self-consistent system of the kinetic equations 
for plasminos and plasmons enables us to have a new look at the known results obtained from the diagrammatic 
technique. It also allows us to represent them in a more visual form and to obtain some new results, in particular, 



1 Here, we have not discussed the problems associated with calculation of the damping rate of moving fermions in a hot QED plasma, 
which possesses proper specific features (see, e.g. Ref. [jlCj]). 

2 Notice that the fact of acquiring the transverse gluons in hot QCD plasmas of the thermal mass was first mentioned by Shuryak |Lq|. 
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connected with the possibility to observe in more detail the pumping dynamics of plasma excitations energy over the 
oscillation spectrum and energy pumping from the fermion branch of excitation into the boson branch and vice versa. 
Our strategy for deriving the relevant kinetic equations is similar to that already used in connection with the purely 
(soft) gauge sector plf . 

From most papers close to the subject of our research the recent work by Niegawa [^3) is to be mentioned. Within 
the framework of the closed-time-path formalism the generalized Boltzmann equation that describes the evolution of 
the number-density functions of fermion quasiparticles was derived by him. The transport equations for "normal and 
abnormal modes" emerge here from the requirement of the absence of large contributions (which is the result of pinch 
singularity) of perturbative scheme proposed in Ref . Q . In spite of the generality of a suggested approach, particular 
expressions for collision terms unfortunately were not given. It makes it difficult to establish the connection with 
other approaches, also dealing with the problem of relaxation of plasma excitations, in particular, with hard thermal 
loop effective theory. 

The paper is organized as follows. In Sec. II, after summarizing the conventions and the notations used in this 
paper, we write a starting coupling set of the soft-field equations and dynamical equations describing the motion of 
hard particles in the presence of the soft quark and gluon fields. The approach to the solution of a given nonlinear 
system based on the approximation scheme method-the weak-field expansion-is proposed. In Sec. Ill, the first-order 
approximation of the induced source is considered, and the correlation functions of fermion and boson excitations 
are introduced. In Sec. IV, the second- and third-order approximations of the induced source and the color current 
induced by the soft-fermion field are studied. In Sec. V, we discuss the consistency with gauge symmetry of the 
approximation scheme used. In Sec. VI we derive the generalized kinetic equation for soft-fermion excitations and 
supplement it with a generalized kinetic equation for soft-boson excitations. The right-hand side of the last equation 
contains the terms responsible for the nonlinear interaction between soft bosonic and fermionic modes. Sec. VII 
presents a detailed consideration of the terms on the right-hand side of the generalized kinetic equation obtained 
in the previous section, their identification with specific physical processes. Then, kinetic equations describing the 
nonlinear interactions between plasminos and plasmons are extracted. In Sec. VIII, by means of the effective Ward 
identities, we show that the function defining the plasmino-hard-particle scattering matrix element squared is the 
same in the covariant gauge as in the temporal one. In Sec. IX, the transformation of the nonlinear Landau damping 
rate for a plasmino similar to the transformation in a pure gauge sector of soft QGP excitations is performed. In 
Sec. X, on the basis of decomposition of the nonlinear Landau damping rate into positive and negative parts, the 
model problem of two interacting infinitely narrow packets with fermion and boson quantum numbers is considered, 
and the kinematic relations between wave vectors of excitations are defined, so that one or another process of pumping 
over of energy occurs. In Sec. XI, an explicit expression for the nonlinear Landau damping rate of the plasmino at 
rest is derived, and a comparison with a similar expression for the damping rate obtained within the framework of 
the HTL approach is carried out. In Sec. XII, the behavior of the nonlinear Landau damping rate for normal-particle 
excitations near the light cone is considered, and this damping rate is shown to diverge owing to the mass singularities 
of quarks. In the Conclusion we briefly discuss a more general statement of the research problem of soft excitation 
dynamics in hot QCD plasma within the framework of kinetic approach outlined in this paper. 



II. BLAIZOT-IANCU EQUATIONS 

We adopt conventions of Blaizot and Iancu p2]| . We use the metric g^ v = diag(l, —1,-1,-1), choose units such 
that c = ks — 1 and note X — (JTo,X), p — (po,p), ■ ■ ■ ■ We consider an SU(7V C ) gauge theory with Nf flavors of 
massless quarks. The color indices for the adjoint representation, a,b, . . . run from 1 to N£ — 1, while those for the 
fundamental representation, i, j, . . . , run from 1 to iV c . The Greek indices a,/3, . . . for the spinor representation run 
from 1 to 4. 

The quark field ip and the gauge field = A^t" 1 with N^ — l Hermitian generators in the fundamental representation 
obey the field equations 

iW,(x) = n(x), (2.i) 

[ly.^^l-rV^W =#(I) 7 ,^(I)t°+j,(I), (2.2) 
where If) = j^D^ is the covariant derivative in the fundamental representation: 

D^d^ + igA^X), 

F^y — F^ v t a is field strength tensor with = d^A^ — d v A<t — g f abc A b A^; [ , ] denotes the commutator and £ is a 
gauge parameter fixing a covariant gauge. 
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The induced source rj{X) on the ri ght- hand side of Eq. (|]l]) can be written as an integral over hard momentum k 
of a density #(k, X) = t a ^A^(k, X) || 

v(X)= 9 J^^m,X), e k = \k\. (2.3) 

A density A"(k, X) is a generalized one-body density matrix mixing fermion and boson degrees of freedom. 

The total induced color current j M on the right-hand side of field equation (^) is a sum of two parts: the current 
j A (X) induced by the soft-gauge field 



it PO = 9 Ij^v^NfiSniiKX) - Sn A (k,X)) + 2N c 6N A (k, X)}, (2.4) 



dk 

(2^) 

and the current j^(X) induced by the soft-quark field 



/Hk 
^v»[6nl(k,X)-6n*(k,X)}, ^ = (l,v), v = k/|k|, 



(2.5) 



where <5n±'*(k, X) = 5n A : <5fa t a are soft fluctuations in the quark and antiquark color densities and 5N A {k, X) = 
SN Aa t a is a soft fluctuation in the gluon color density. The superscripts indicate the nature of the background field 
that induces fluctuations. On the space-time scale (gT) -1 these functions and density $.(k, X) satisfy the following 
system of equations: 

[v-D x ,6n A : (KX)]=T9vE(X)^&, (2.6) 

ae k 

[v ■ D Xl SN A (k, X)} = — gv ■ E(X) ^r^, (2.7) 

de k 

[v ■ D x ,8nt(k, X)} = ±^- t a ({p(X)t a M(k, X) - £(k, X)t a ^(X)), (2.8) 

(vD x )A(k,X) = -tgC F [N(e k )+n(e k )}^(X). (2.9) 

Here, E(X) = F, a (X)t a is the chromoelectric field, E i = F l0 ;n(e k ) = l/(exp(e fe /T) + l) andiV(e fc ) = l/(exp(e fc /T)-l) 
are the fermion and boson occupation factors, T is the temperature of the plasma, and Cf is the quadratic Casimir 
invariant of the fundamental representation. The function $(k, X) is related to $.(k, X) by $(k, X) =feJ(k, X)j°, 
where the Hcrmitian conjugation refe rs to c olor and spinor indices. 



The self-consistent system of Eqs. (2.1) - (p.9|) for the soft fluctuating fields ip and and their induced sources 



was first derived in Ref. |22j. In the subsequent di scus sion it will be called the Blaizot-Iancu equations. However, as 



distinct from the original paper [[22|, we take Eqs. (2.6) - ( |2.9| ) not as kinetic equations, i.e., time-irreversible ones, but 
as exact "microscopic" dynamical equations coming from the hard thermal loop effective action and describing the 
evolution of the collisionless plasma with zero expectation values of the soft fields (or the associated induced sources). 

The Blaizot-Iancu equations are solved by the approximation scheme method - the weak-field expansion. For this 
purpose first we expand soft fluctuations of the (anti)quark and gluon color densities as power series in the oscillation 
amplitudes of the functions and ip, 

oo oo 

5n A = J2Sni (n \ 8N A = Y. SNA(n) > ( 2 - 10 ) 



5nt = y j 5nt^ 2 \ (2.11) 



7X = ^vX (n ' 1) , (2.12) 
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where the index n shows that 8ri±^ n \ 8N A ^ n \ 5n^} n ^ and A^™' 1 ) are proportional to the nth power of A^. Since 
the fermion fields appear in explicit form on the right-hand side of Eqs. ( |2.S|) and ( |2.9| ), the functions 5n*^™' 2 ^ and 
^(n,i) & \ so depend on amplitudes ip and ip. By virtue of the structure of the right-hand side of Eqs. (2.8) and 
we have 



6n 



*(n,2) 



n,l) 



for arbitrary values of n. 



The induced color currents (2.4) and (2.5) and the induced source (2.3) are expressed as 

oo 



(2.13) 



where 



and 



J/J, / , J/J, 



n=0 



dk 



v^[N f {5n^ n) - 5n A(n) ) + 2N C 8N A ^}, 



dk 



(27T) 



*(n,2) c *(n,2) n 

' — on 



(2.14) 

(2.15) 
(2.16) 



where 



" = E ? ? (n,1) ' 

n=0 



(2tt)3 e k 



(2.17) 



(2.18) 



Now we turn to the field equations (2.1) and (2.2), connecting the quark and gluon soft fields with induc ed s ource 
r] and total color current j^. Let us rewrite these equations, explicitly separating the free parts in Eqs. (2.1) and 
( 2T2I) from interaction terms. Taking into account the expansion of induced current ( 2.13| ), ( |2.14 ) and induced source 
( 2 . 1 7[) , we have 



(0,1) _ 



gfiip + rj NL , 



(2.19) 



f NL +g$l v t a ^t a -igd»[A^A v ]- ig{A^(F^)L}+ g^A^A^A"}}. (2.20) 



Here, the indices L and NL denote the linear and nonlinear parts of strength tensor, the induced source and the 
color-induc ed cu rr ent w ith re spect to A^ and tp. The approximation that will be made is the truncation of the series 
expansion ( l2T3| ), (pT4j) and J2T7| ). The accuracy of this approximation is controlled by the characteristic amplitudes 
of the soft fields, which will be discussed in Sec. V. As will be shown to account for the nonlinear interaction between 
waves and hard particles in QGP at leadin g ord er in q, it is su fficie nt to restrict the consideration to the third order 
in total powers of ip and A^ in expansions (2.13), (2.14), and (2.17). 



III. LINEAR APPROXIMATION OF THE INDUCED SOURCE n. CORRELATION FUNCTIONS OF 

THE FERMION AND BOSON EXCITATIONS 



We now come to the derivation of the kinetic equation for the soft- fermion modes. The starting point is Eq. (2.19) 



Its left-ha nd side contains the linear approximation of the induced source 77, whose explicit form can be easily found 



from Eq. (2.9). We prefer to work in momentum space. The corresponding equations are obtained by using 
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and similar transformations for Sn±, SN A , etc. Here and in what follows we denote the momenta of the soft-quark 
fields by q, q', qi, . . . a nd th e momenta of the soft-gauge fields by p, p' , pi, .... 

Let us linearize Eq. (2.9) with respect to ip and A^. The result of the Fourier transformation of the linearized 
equation is 



' 1 \k,q)=gC F 



v ■ q + it 



[N{e k ) + n{e k % e +0. 



(3.1) 



By substituting Eq. (3.1) into the Fourier transform of relation ( 2.18 ) (for n = 0) and performing the radial integration 
over ciefc, we obtain the linear approximation with respect to the fermion field of the induced source, 



r 1 ^ 1 \q)=SE(q)^q), 



(3.2) 



where 



q + ie 

is a well-known HTL expression for the soft-quark (retarded) self-energy, Uq = g 2 C F T 2 /8 is the plasma frequency of 
the quark sector of plasma exci tation s, and d£l is an angular measure. 

Furthermore, we rewrite Eq. (2.19) in momentum space. Taking into account Eq. (|3.2|) we obtain 



*S-\q'W{q') = - g U(piW(qi)6(q' - q 1 -p 1 )dq 1 dp 1 - V) - V ^ j (q'). 



(3.3) 

Here *S(q) = [— 4+ S'Ejg)] -1 is an equilibrium propagator for a soft quark, corrected to leading order in g. 

Let us multiply Eq. ( |3.3| ) by the Dirac conjugate amplitude ^(—q) = (ip l (q))^° and take an expectation value over 
the off-equilibrium ensemble: 



m-q)*S-\q'W(q')) =-g (r(~q)(MpiWqi)) J ) S(q' - qi - Pi)d qi d Pl 



(3.4) 



- m-q)v (1 ' 1)j (q')) - m-<l)v i2 ' 1)j (q'))- 
We now introduce the correlation functions of soft-fermion and -boson excitations 

nu?'?) = m-iw<i')}, i%(p',p) - (a^waKp)), (3.5) 

respectively. The asterisk denotes a complex conjugate. The considered soft excitations are necessarily colorless and 
have zero fermion number by virtue of the fact that mean fields (AV) and (ip l a } or the associated mean induced color 
current and source a re a ssumed to be vanishing. Therefore, for the physical situation of interest, the off-equilibrium 
two-point functions ( |3.5| ) are diagonal in color space, which will be implied in what follows. 

For the conditions of stationary and homogeneous QGP's (i.e., when correlation functions ( |3.5D in the coordinate 
representation depend only on the difference of coordinates and time AX = X' — X), we have 

The QGP state becomes slightly heterogeneous and nonstationary because of the effects of the nonlinear interaction 
between waves and particles. This leads to a (5-function broadening, and T 3 g a and J~J depend on both arguments. 

Let us introduce T^(g', q) = T* a (q, Aq), Aq = q' - q with | Aq/q \ « 1 and lf v {p' ,p) = I${p, Ap), Ap = p' - p 
with | Ap/p | <Cl, and insert the correlation functions in the Wigner form 

T&(?.*) = / Tg,(g, AgJe-'^-'dA?, T*{j>,x) = f I$(p,Ap) e- i ^ x dAp, (3.6) 



slowly depending on x. In Eq. (3.4), we replace q <-> q' , i <-> j, take a complex conjugation, and then subtract the 
resulting equation from Eq. (3.4), expanding beforehand the quark self-energy into "Hermitian" and "anti-Hermitian" 
parts: 

*E(g) = 6^(q) + 5Z A (q), 7 (<ffi H (g)) V = SZ^q), 7 (<5S A ( g ))t 7 = - 5Z A (q). 



7 



We assume that the anti-Hermitian part ST, A is small relative to the Hermitian part 5S H and is of the same order as 
the nonlinear terms on the righ t-hand side. We can therefore set 5T, A (q) ~ S'E A (q') and move the term with ST, A into 
the right-hand side of Eq. (3.4). We expand the remaining terms in the left-hand side with respect to Aq up to the 
first order. This corresponds to a gradient expansion procedure usually used in the derivation of kinetic equations. 
Multiplying the resulting equation by e- tA i- x an( } integrating over Aq, we obtain 



tr (^ [ -^ + ^ H(9)] 



dT^(q,x) 
dx u 



2*tr (Sll A (q)T ji (q,x)) 



(3.7) 



-igjdq'dq 1 dp l {{^{-q 1 )^{-p 1 )) i ^{q'))5{q-q 1 -p 1 ) - (V7(-g)(4(Pi)^i)) J )%' - Qi ~ Pi)} 

dq'{((n {1 ' 1)l (~qWW)) ~ $*(-9)»7 (1,1)i (9')» + ((v {2A)l (-q)^(q')) ~ (r(-q)v (2A)j (q'))j\- 

Here, the Dirac trace is represented by tr, fj(—q) = f^(<z)7°, and we take into account the reality of the gauge field: 
A*^ a {p) — A®(—p). The linear term with <5£ A on the right-hand side corresponds to the linear Landau damping of 
soft-fermion excitations. 



IV. SECOND AND THIRD APPROXIMATIONS OF THE INDUCED SOURCE rj 



Now we are concerned with computation of the nonlinear corrections to the induced source on the right-hand side 
of Eq. (3.7). At first, we define v)^ 1,1 '. Substituting the series expansion (2.12) into ( |2.9| ) and keeping only the terms 
of the second order in fields, we find 

(v ■ d x )^(k,X) = -ig(v ■ A(X))^(k,X). 

Performing the Fourier transformation of the last equation and taking into account an explicit form of ^°' 1 ^ > ( |3.l| ), 
we derive 



A^^k, q) = g 2 C F [N(e k ) + n(e k )} - t t a 



1 



v ■ q + it J v ■ qi + ie 



■ A°(pi)ilj(qi)S(q -qi- pi)dqidp x 



(4.1) 



Substituting the obtained expression into Eq. ( 2.1q ) (for n = 1) and performing the radial integration over de k , we 
find the required induced source correction 



r,^(q)=g^t"l^ 



an 



47r (v ■ q + ie)(v ■ q\ + ie) 



AZip^ipiqx) s (q -qi- p^dqidpx 



(4.2) 



The expression for the induced source of the third order in the fields is defined by means of reasoning similar to the 
previous one. Here, we have 

T 7 ( V~~T( Im ^(Pi)4(PaM«Zi)<f(3 - ?i - Pi - P2)dq 1 dp 1 dp 2 . 

Ait {v ■ q + ie)(v ■ (q - pi) + ie)(v ■ qi + ie) F 

(4.3) 



Now we return to the initial equation for the soft-fermion field (3.7). We substitute nonlinear corrections to the 
induced source p (E qs. (4.2) and ( |4.3| )) into the right-hand side of this equation. After simple algebraic transformations, 
instead of Eq. (3.7), we find 



U-[-^[-d+SX*(q)} dr 2 q,X) )=2/.r.^i< / )T"(<i-.-,) 



i 5 /"^^{(^(-g)^ 



m-qi)*T^ a (-p 1 ;-q 1 ,q)yA*^(p 1 )^(q'))S(q- qi - Pl ) 



(4.4) 
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Here. 



+ ig 2 j dq'd qi d Pl dp 2 { (A a » ( Pl )A b » {p 2 W (~q) (ST<$> ab (- Pl , ~p 2 ;q' ,~qiMqi)y)S(q' - Ql - Pl - p 2 ) 

T«> a (p; qi ,q 2 ) = t a T<W(p; qil q 2 ) = t a (^ + ST™ (p; qi ,q 2 )) (4.5) 

is an effective (i.e., HTL-resummed) vertex between a quark pair and a gluon, that is, a sum of the bare vertex 7 M 
and a corresponding HTL correction J|, [2^, [24|] 

«T«Wfc) = -, (4.6) 

J 47r (v ■ qi + ie){V ■ q 2 — ie) 

and 

<* r ^ (Pi.Pa;gi,?2) = -w / — — : — — , — - — — — + — -, — - — r— - (4.7) 

J 47T (u ■ 91 + ie)(v ■ q 2 — ie) \v ■ (qi + pi) + ie v ■ (qi + p 2 ) + ie J 

is an effective vertex between a quark pair and two gluons (this vertex does not exist at tree levei, and in the 
leading order it arises entire ly f rom the HTL §, H|). The superscript (Q) denotes that the vertex ST^X; Yi, Y 2 ) 



corresponds to the function (4.6) in the co ordin ate representation, where the time arguments satisfy > X Q > Y 2 ° 



(boundary conditions), the vertex function (4.7) corresponds to ST^ (X\ , X 2 ; Y\,Y 2 ), whe re the time arguments satisfy 
Yi > X® > X® > Y§ for the first term in parentheses on the right-hand side of Eq. Q4.7[ ) and for the second term we 
have YP > X® > X® > Y§, i.e., the time argument of the external quark leg incoming in the vertex functions (4.6) 
and (p~7|), is largest. In deriving Eq. (4.4) we have dropped the terms proportional to (A^ip) (A u ip) not contributing 
to the right-hand side of kinetic equation for the soft-fermion excitations at the leading order in g, in which we are 
interested. 

Because of nonlinear wave interactions, phase correlation effects occur. By virtue of their smallness, the expectation 
value of a product of four quantities {A ll A v tp^i) can be expressed approximately as a product of the expectation values 
of two fields (A^A U ) and (^ip). For a product of three fields this approach yields zero; in this case, the weak correlation 



between the fields is to be taken into account. The third-order correlation functions on the right-hand side of Eq. (4.4) 
contain amplitudes of waves with different statistic, nonlinear corrections to that are defined by corresponding field 
equations. Let us consider first of all the nonlinear correction to a free quark field. For this purpose, we use field 
equation keeping in the right-hand side only the second-order terms (with respect to A^ plus iff) 

*S- 1 (q)i/j(q) =-gj T») "(p i; qi, -gM„(piM«l) S(q - gi - Pi)d qi d Pl . 
The approximate solution of this equation has the form 

i>{q) = ^ ] {q) - 9*S(q) J TW)"( Pl ; 9l , -q)A^ ( Pl )^°\ qi ) 5(q -q x - Pl )d qi d Pl , (4.8) 
where A^ and ip^ are solutions of the appropriate homogeneous field equations corresponding to free fields. 



Now we consider the nonlinear correction to a free gauge field. For this purpose, we use field equation (2.20) is 
rewritten in momentum space, keeping on the right-hand side the terms of the second order in fields: 

•V- 1 ^(p)AZ(p)=j A W**(p)+j*lW a *(p)+g[ {j(- qi )rt a i,( q2 )S(p + qi - q 2 )d qi dq 2 



- ^9f abc lr^\p,- Pl ,- P2 )Ai( Pl )Al( P2 )5( P - Pl - P2 )d Pl dp 2 . (4.9) 

Here, *T> tLV {p) — — [p 2 g tlL/ — (1 + £~ 1 )p >i p v + SH fl ^(p)]~ 1 represents the medium-modified (retarded) gluon propagator 
with the soft- gluon polarization tensor, 

1 \ J 47T v ■ p + ieJ 
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corrected at the leading order in g, oj^ — g 2 T 2 (Nf + 2iV c )/18 is a plasma frequency of the gauge sector of plasma 

exci tatio ns, and T flvX (p,pi,p 2 ) is a bare three-gluon vertex. The first and the last terms on the right-hand side of 
Eq. ( |4.9| ) give the contribution to a correlation function {A'tptp) proportional to (AA'tpip) — (AA) (V>^/>) . By replacing 
interacting fields by free ones and dividing the fourth-order correlation into a product of the second-order correlators, 
this contribution vanishes. Therefore these terms in this approximation c an be dropped. 

Now we consider the remaining terms on the right-hand side of Eq. (4.E). For determining the explicit form of 
^.*(o,2)o ^ j g necegsar y ^ derive the expression 5n^S°' 2 " >a according to Eq. ( 2.16| ). For this purpose, we substitute the 
expansion (2.11) into Eq. ( |2.8| ), keeping only the ter ms w ith quark fields V>, ip on the right-hand side. Performing the 
Fourier transformation and taking into account Eq. (3T), we find 



*(0,2) a. 



1 



l (k,p) = ±g 2 C F — [N(e k ) + n(e k 



1 



(v ■ qi — ie)(v ■ q 2 + ie) 



— ip{-qi) ^t a tp(q 2 )S(p + qi - q 2 )dq\dq 2 . 



Substituting the obtained expression into Eq. ( 2.16| ) and performing the radial integration over de k , we find the 



required current correction j*^ 0,2 ' : 



J 



*(0,2)a/ X _ 



(p) = - .9' 



ip(-qi)jtt a ip(q 2 )S(p + qi- q 2 )dqidq 2 



An (v ■ p + ie) \v ■ q 2 + ie v ■ q± — ie / 
Taking into account the discussion above, we obtain the following equation, instead of Eq. (|4.9|): 



(4.10) 



*V- 1 ^(p)A a u (p) = g i,(~q i rT^(p;q ll ~q 2 )t a 1 p(q 2 )S(p + q 1 - q 2 )d qi dq 2 . 



(4.11) 



Here. 



(4.12) 



dn 



5Tf ] {p]qi,q 2 ) = -u% 

J 47T [v ■ qi - ie)(v ■ q 2 - le) 

is the effective two-quarks-one-gluon vertex function, where now (in the coordinate representation) the time arguments 
satisfy X° > max (Yf, Y 2 °) and the chronological order of Yi and Y 2 is arbitrary (2^] . The time argument of an ex terna l 
gluon leg incoming in the vertex is largest, as indicated by superscript (G). The approximate solution of Eq. (4.11) 
(with accuracy required for our further calculations) is of the form 



A;(p) = 4 0) °(p) + 9 *2V(p) / ^ (0) (-9l) *^ G ^( P] qi , -q 2 ) t a ^°\q 2 ) 5{p + q 2 )d qi dq 2 . (4.13) 



At the end of this section we present the expression for the next term in the expansion of current j*, which is 
needed for deriving general kinetic equation for soft-gluon modes in Sec. VI. Performing similar calculations, we obtain 



f, (1 ' 2)a (P) 



2 a fdQ v^(f) aP 

9 w o / T~ ~~ — 

' 47r v ■ p + ie 



(T a ) bc {t c y 



(v ■ qi - ie){v ■ q 2 + ie) 



a,b\ij _ 



(t a t") 



1 



(v ■ {q 2 +pi) + ie)(v ■ q 2 + ie) 



- (t b t a ) tj 



1 



(v ■ (qi - pi) - ie)(v ■ qi - ie) 



(4.14) 



* A l{PiW a (-<liWB{<l2)5{p- Pi +qi - q 2 )dp 1 dq 1 dq 2 



where (T a ) bc = —if ahc . The appearance of the term with the production of generators of the SU(A r c ) group in other 
representations (the adjoint and the fundamental ones) is a special feature of the last expression. This term is not 
vanishing only in the processes of the highcrordcr in field powers, the processes of Boltzmann type, i.e., the scattering 
processes between the soft-fermion and -gluon excitations (see Conclusion). 
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V. CHARACTERISTIC AMPLITUDES OF THE SOFT FIELDS 

In this section we will estimate the typical ampli tudes o f the soft fie lds, b oth in coordinate and in momentum space, 
such that the truncation of the series expansions ( [2.13 ), ( [2.14 ), and ( [2.17 ) can be made and in the long run one can 
derive a closed system of gauge-invariant kinetic equations for soft plasma excitations. 

For this purpose, let us discuss in more detail the approximation scheme used in this paper. In fac t, t he value 
6 = gRA^ is a (dimensionless) parameter of the expansion in powers of nonlinearities in Eqs. (2.1) - (|2j|), where 



R ~ d x is a typical length for off-equilibrium deviations. It follows, for example, from rewriting the equations of 



motion for the soft fields Q2.1| ) and (2.2) in the form (2.19) and (2.20), where the gauge fields A^{X) are expanded 
from the covariant derivatives — + igA^{X) and the fiel d stre n gth te nsor F,. V (X ). The disadvantage of the 
proposed iterated scheme for deriving the terms in expansions ( 2.13 ), ( 2.14 ), and ( |2.17|) is exp licit breaking of the 
non-Abelian gauge symmetry of the dynamical equations for the soft fluctuations fl2.6| ) — (2.S) at each step of the 
approximate calculations. Here, the answer to the nontrivial question on the recovery of a gauge symmetry depends 
directly on the order of the magnitudes of the soft fields, the dynamics of which we would like to describe. Below, we 
discuss two relevant cases. 

Let us assume that the amplitudes of soft-boson and -fermion excitations are of orders g coinciding with the orders 
of An and %jj presented in Ref. E3], i.e., 



\A»{X)\~T {\A»{p)\~l/g{gTf), 



mx)Mx)\~gT* (\mU(q)\ ~i/g\gTf). 



In this case, by using the obtained expressions for the terms in the expansion of the induced source r\ ( |3.2| ), (4.2), and 
(4.3) and the color currents |H| and j* (4.1C) and (4-14), we have the following estimates: 



^%)~r/M) (g) ^(2,i) (9) 



gT |V(«)|, 



1 

7r' 



g(gT) 4 $(q)Mq)\ 



g 2 T 



Thus, all terms in each expansion ( 2.13| ), (2.14), and (2.17) are of the same order of magnitude and the problem 
of resummation of all the relevant contributions arises. Thus, gauge symmetry is recovered. In a pure-gauge case 
such an approach is represented, for instance, in the work of Jackiw and Nair |2q] dealing with the derivation of the 
non-Abelian version of the Kubo formula. Here, the explicit expression for the induced current j^(p), including the 



contribution of all higher point functions, that actually coincide with the iterations j^ y ' -', ... |21|, is defined. 
The requirement of gauge invariance leads to the completely nonlinear theory and principal impossibility of deriving 
closed kinetic equations for number densities (i.e. two-point functions) of the soft-fermion and -gluon excitations. 

In this work, we restrict consideration only to a finite number of terms in expansions of the induced source and 
color currents. This imposes more rigid restrictions on the order of the magnitudes of soft fields. Blaizot and Iancu 
p6| (Sec. 1.1) showed in the special case when the soft fields were thermal fluctuations at the soft scale gT (namely, 
this situation takes place in our case), their typical amplitudes would be of the orders 



IWI-fl^T (\A,(p)\ 



1 



9 1/2 (gT) 3 



mx)\mx)\~g A T* m q )\m q )\ 



g{gT)\ 

In this case the parameter of expansions 6 ~ g 1 ^ 2 <C 1 and we have the estimates 

V W(q) ~ gT\^(q)\, n^(q) ~ <? 3/2 T V ^(q) ~ g 5/2 T\^(q)\, 



11 



;A (2) 



^ (0 ' 2) b)-3(5T) 4 |^(<7)||^(g)| 



. i* (1 ' 2) (p)~5 3/2 (5T) 4 |Vi(g)||V(g)| 



1/2 T 



i.e., every successive term in the source and current expansions is suppressed by one power of g 1 / 2 in comparison with 
the preceding ones, and the use of the perturbation theory is therefore justified. By taking into account the nonlinear 
interaction betwee n soft fields and hard particles at leading order in g, i t is s ufficient to keep only the first three terms 
in the expansions (2.13) and (2.17), and first two terms in expansion (2.14), as mentioned at the end of Sec. II. The 
recovery of a gauge symmetry here, occurs another way, by the account of weak correlation in the calculation of the 
expectation value of three soft fields in a slightly heterogeneous, slightly nonstationary quark-gluon plasma (see the 
following section). In this case we derive a closed system of the gauge-invariant kinetic equations that describe the 
evolution of the off-equilibrium deviations in the number densities, and obtain the damping rates from the collision 
terms, which are closely allied in form to the corresponding damping rates in the HTL approximation S, 0]. 



VI. GENERALIZED KINETIC EQUATION FOR SOFT FERMION EXCITATIONS 



Now we turn to the equation for the soft-fermion excitations (4.4). We substitute Eqs. (4.8) and ( [1.13 ) into 
third-order correlation functions entering into Eq. (4.4). Because of the fact that "4>(°\ ip(°\ and Affl represent 
the amplitudes of entirely uncorrelated waves, the correlation function {A^ ipWi/jW) drops out. In this case, each 



term in (i} l {-q)A a ^ l {p 1 )( *T { ^> a ^{qi)) j ) and {{^{-qi)*Yf )a ) i A* a ^{p 1 )^>{q')) should be_ defined more ex_actly. In 
the fourth-order correlation functions, within the accepted accuracy, we replace the fields i/j, i/j, and A^ by i/j^ ', 

and A { °\ 

Furthermore, we make the correlation decoupling of the fourth-order correlation functions in terms of the pair 
correlation functions by the rules 

$W (-qi)^ 0) (g 2 )^ (0) (-<Z 3 )V> (0) (<Z4)} - (V^ (0) (-<Zl)V> (0) (<&)) (^ (0) (-Qs)^ (q 4 )> 



',(0) 



(0), 



'.(°>f 



(A^(p 1 )A^(p 2 )^ Q \-q 1 )^°\q 2 )) = (A(°)(p 1 )A( )fc))(^°)(- ( ji)^ 0) (g2)). 

Here, we suppress the color, spinor, and Lorentz indices, and take into consideration the Grassmanian nature 
of the quark field. Setting T^ Q = P'T^ and = 5 ab I^ Ul taking into account 5 ST^~,(jpi,p2;qi,q2) = 
Cf&Y nv{pi,P2'i qi, 92), after cumbersome calculations, we come to the following generalized kinetic equation 3 for 
the soft-fermionic excitation QGP, instead of Eq.( |4.4D 



tr 



_d_ 



m n {q)] ^^Y\ =2 itr(5Z A (q)r( q ,x)) 



(6.1) 



+ 2g 2 C F [ d qi lm 



*V, u (q - gi ) tr{x(g) *r«>"(g - q x ; q u -g)T(g x ) *T^^(q - q x ; q 1: -g)} 



d gi d Pl V(pi) Im[tr{*rW)^(-p i; g)*5(g)*rW) "(p i; g lf -g)T(gi)}] S(q -q x - Pl ) 



3 The term generalized here means that we have not yet restricted ourselves to any mass-shell conditions. 
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-Jdpl^p) Im[tr{ (jrW^Cp, -p; g, -g) - T«> "(p; g - p, -g)*S(g - p)*r«>"(-p; q, -q+p) 



T«> "(-p; q + p, -q) *S(q + p) T«> "(p; g, -g - p))x(g) } 



On the right-hand side of Eq. ( p.l| ) the x dependence of Yp a (q) and I^vip) is understood, although not explicitly 
written. 

In deriving Eq. (^^) we use the following properties of effective tw o- q uar k- one - gluon and two-quark-two-gluon 
vertex functions, which immediately follows from initial definitions ( |4.5| ), ( [4.6| ), (4.7) and ( 4.12| ) 



7 *r^)t(p; gi ,g 2 ) 7 °= *rW)(p;g 2 ,g 1 )= T«)(-p; -g 2 ), 
7° Tpfe 91,^)7°= *rf)(-p;-g 1) -g 2 )= *Tf\-p;-q 2 - qi ), 

1 8Tf l )\p llP2 -q 1 , q2 ) 1 a = - <5rW>(pi, pa; fc, gi ) = - 5rW)(-pi, -pa; 

Moreover, we assume that under Hermitian conjugation the Wigner function T(q, a;) behaves like an ordinary 7 matrix: 

7°rt(g,x)7° =T(q,x). (6.2) 

Below this p rope rty is shown to ensure that the physical variable such as the fermion number density is real. 

Equation (6.1) is incomplete, since the unk nown function I tlv {p 1 x) en ters this equation. This function obeys the 
kinetic equation, which is similar to Eq. (6.1). The soft-field equation ( 2.20 ) is initi al for deriv ation of this kinetic 



equation. We use the obtained expressions for the induced color current j* ( 1.10 ) and ( 4.14 ) and corresponding- 



expressions for the induced color current defined in Ref. |2lJ . Performing calculations similar to previous ones and 
keeping in the right-hand side only the terms responsible for the nonlinear interaction between soft-fermion excitations 
and soft-gluon ones (purely gauge sector was considered in detail in Ref. pi|), we come to the following generalized 
kinetic equation for soft-gluon excitation QGP's (again, the x dependence of most quantities is suppressed) 



_ (i + r i yy + su^{p)] dI ^ x) = - 2i <m A ^(p)V(p, x) 

- 2g 2 T F (lm( *V^(p)) J d qi dq 2 tr{T(g x ) "(p; Ql , -g 2 )T(g 2 ) *r( G ) "(-p; -q x , g 2 )} 
+ V(p) /dgIm[tr{(5r( G ) ^(p, -p; g, -g) + T<<« "(- PlP - q, q )*S(p - g)*r( G ) "(p; -p + g ,-g) 



(6.3) 



(p; g, -p - g) *S(p + g) T«> "(-p;p + g, -g)) T(g)} 



Here, Tp is the index of the fundamental representation, 5H HlJ,v (p) and £n A ' tly (p) are the Hermitian and anti- 
Hermitian parts of the soft-gluon self-energy, respectively, and ST^ M " is the effective two-quark-two-gluon vertex 
function defined by an expression 



6T^^( Pl , P2;qi , q2 )=-^f^ — 

J 47T (V • 1 



1 —, \ H 7 ; , (6.4) 

gi - ie)(v ■q 2 -ie)\v- (qt+p^ + ie v- (gi +p 2 J - ze/ 

where the first term in large parentheses satisfies the chronological order: X® > Y 2 , X® > X 2 > Y® in the coordinate 
representation, and the second term satisfies the chronological order X° > Yj , X® > X® > Y 2 \ i.e., the time argument 
of one of external gluon legs coming into the vertex is largest. The mean of the superscript (G) is ev ident here. Notice 
that in deriving Eq. (|6.3|) we use the properties of Hermitian conjugate of the vertex function (6.4): 



7 (5r^ )t (p 1 ,p2;<?i,g 2 )7 = ST ( j ^>(p 2 ,p 1 ; qi , q2 ) = - $Tffi(p u p 2 ; q 2 , qi). 

In closing this section we again call attention to t he chrono logi cal order of arguments in the vertex functions entering 
into the right-hand side of generalized equations ( |3.l| ) and ( |6.3[ ), which appear here in nonevident manner. Also, we 
call attention to the fact that the righ t-hand side of Eq. ( |6.5| ) is proportional to the third group invariant Tp as distinct 
from the right-hand side of Eq. (S.l), proportional to the quark Casimir invariant Cp and the kinetic equation for 
l^y (p, x) with purely gauge sector on the right-hand side pl[ , proportional to the gluon Casimir invariant Ca ■ 
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VII. CLASSIFICATION OF TYPES OF NONLINEAR WAVE INTERACTIONS. PLASMINO AND 

PLASMON KINETIC EQUATIONS 



Now we perform a preliminary analysis of the right-hand sides of the generalized kinetic equations ( |6.l| ) and ( p.q) , 
the purpose of which is the identification of the terms with specific physical processes described by them. To establish 
this connection it is convenient to use the Feynman diagrams, defining the contribution of the leading order in the 
coupling constant to the damping of soft-fermion excitations in framework of resummed perturbation theory. There 
are known || to be two diagrams (Figs. 1(a) and 1(b)) with soft loop momentum that contribute at order g to the 





(a) (b) 

FIG. 1: One-loop diagrams that contribute to the imaginary part of the self-energy of quarks. The blob stands for HTL 
resummation. 

effective self-energy for quarks: the usual self-energy graph at one-loop order (Fig. 1(a)), except that all of the vertices 
and propagators are effective, and the effective gluon tadpole on the soft-fermion line (Fig. 1(b)), which is special to 
the effective expansion. 



Let us consider the second term 4 on the right-hand side of Eq. (6.1). Physically, this contribution just corresponds 
to the stimulated scattering processes of soft-fermion excitations by hard-particle QGP through a soft-virtual gluon 
oscillation, without the change of statistic type of both soft and hard excitations, i.e., 

q + Q(Q)^q + Q(Q), q+G^q + G, 

and similarly for q. Here, q represents fcrmionic collective excitations (we do not distinguish between normal-particle 
excitations and plasminos) and Q, Q, and G are quark, antiquark, and gluon excitations with characteristic momenta 
of order T. It can be easily inferred if cutting the usual self-energy graph before spelling out the effective gluon 
propagator, i.e., inserting a hard bubble along the gluon line, as depicted in Fig. 2. 



The third term on the right-hand side of Eq. (6.1) is associated with the so-called decay processes depicted in Fig. 3 



following from the cutting of the usual self-energy graph 

q q + g, q + q -> g, 

where g is the gluonic collective excitations (here, we do not distinguish between longitudinal and transverse excita- 
tions). Similar decay processes have been studied in detail, since they are immediately connected with the processes 
of soft-dilepton productions in hot QCD plasma. In the high-temperature QCD plasma in the approach based on 
using dispersion relation methods [p7[, these processes are connected with the so-called pole-pole terms. 



Finally, the remaining terms grouped on the right-hand side of Eq. (6.1) are connected with the more interesting 
process of the nonlinear wave interaction with stimulated scattering of soft-fermion excitations by hard-particle QGP's 
varying with the change of statistics of excitations: 

q(q) + G^g + Q(Q), q(q) + Q(Q)^g + G. 

It is defined by two different physical processes. The first represents the Compton scattering type process and is 
connected with the eff ectiv e two-quark-two-gluon vertex function (more exactly, with its imaginary part) on the 



right-hand side of Eq. (OA). This corresponds diagrammatically to the appropriate cutting of the gluon "tadpole" 
graph, before illustrating the effective two-quark-two-gluon vertex function, as depicted in Fig. 4. The second process 
connected with the remaining terms defines the scattering of a fcrmionic quantum oscillation through a virtual 



4 As was mentioned at the end of Sec. Ill, the first term on the right-hand side of Eq. (6.1), involving the imaginary part of the HTL 
self-energy, describes linear Landau damping. It vanishes, when one studies the decay of on-shell excitations and therefore can be 
omitted. 
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FIG. 2: The process of the stimulated scattering of soft fermionic excitations by hard QGP particles through a resummed gluon 



f G) 

propagator T>, where a vertex of a three-soft-wave interaction is induced by 1^ . The double lines denote hard particles 





(b) 




FIG. 3: (a), (b) Radiation (absorption) of the soft-boson excitations by soft-fermion ones; (c) annihilation of the soft-fermion- 
antifermion excitations into soft-boson excitations. 



soft-fermion oscillation. Diagrammatically this process of scattering is denned by the corresponding cutting of the 
usual self-energy graph before illustrating the effective quark propagator as in Fig. 5. As it will be shown in the 
next section, the terms corresponding to the processes drawn in Figs. 4 and 5 form a gauge-invariant function. The 
process, described by these terms, will be called further the nonlinear Landau damping of the soft-fermion excitations 
by virtue of the great similarity of the expression defining this process and the corresponding expression in a purely 
gauge case Ml- 



Here, it is necessary to note that the terms on the right-hand side of kinetic equation ( 3T ) corresponding to the 
process diagrammatically shown in Fig. 5 are very close in structure to the terms connected with stimulated scattering 
processes, without changing the statistic type and which are diagrammatically given in Fig. 2. However, in the last 
case the Compton scattering type contribution is absent, which, in turn, is associated with the absence of the effective 
quark tadpole on the soft-quark line. For this reason, we not assign this process to the one of the nonlinear Landau 
damping. 

Now we turn to the right-hand side of generalized kinetic equation J6.3| ). The second term on the right-hand side 
of Eq. (6.3) is defined by time reversal of the decay processes that are drawn in Figs. 3(a) and 3(c). The remaining 
terms on the right-hand side of Eq. (6.3) are associated with the processes defined by the time reversal of the processes 
depicted in Figs. 4 and 5, i.e., with the process of the nonlinear Landau damping of the soft-gluon excitations. 
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FIG. 4: The Compton-like scattering of soft-fermion excitations by a hard QGP particle with a change in statistics of soft 
excitations (s channel). 



I 
I 
I 




FIG. 5: The process of the stimulated scattering process of the soft-fermion excitations by hard QGP particles through a 
resummed quark propagator *S, where a vertex of a three-soft- wave interaction is induced by (t channel). 

In the remainder of this work, we restrict our attention to the detailed study of the process of nonlinear Landau 
damping, and moreover we restrict ourselves only to the process of the nonlinear interaction between plasminos and 
plasmons, i.e., purely collective excitations in hot QCD plasmas. The exception is only Sec. XII, where we consider 
the difficulties connected with the calculation of the nonlinear Landau damping rate for normal quark excitations in 
the vicinity of the light cone. 

As it was mentioned in the Introduction in global equilibrium QGP's for the case of massless quarks and for zero 
chemical potential there are two branches of fermion excitations with a positive energy. In this connection we define 
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the Wigner function T(g, x) in the form of an expansion 

T{q,x) = h+(q)f + (q,x) + fe_(q)T_(g,x), 



(7.1) 



where T± are certain scalar functions and h± (q) = (7 =F q • 7)/2 with q = q/ 1 q| . By the condition ( 0/2 ) the functions 
T± are real. Let us define dependence on variable q°. 

For this purpose we omit nonlinear terms and the anti-Hermitian part of the quark self-energy in Eq. (3.4). Substi- 
tuting the function S : > l (h±(q})p a T±(q, x) S(q — q') instead of (— q^^q')} , we find 



where 



with 



-[ q° T (|q| + Re 5E±(q))]T± (q,x)=0 
8Z ± ( q ) = ±~tr{h±(ci)5X( q )} = ^ 



(7.2) 



In 



l + z 



(i-N) 



The solutions of these equations have the structure 

T±(<7, x) = T±(q, x) S(q°- w±(q)) + T T (-q, x) <%° + w T (q)). 



(7.3) 



Here, T±(q, x) = T q are certain functions of a wave vector q and uj±(q) = w q are frequencies of the normal-particle 
excitations and plasminos, respectively. The solutions (7^3) describe what is called a quasiparticle approximation. 
In this approximation we assume that the off-equilibrium two-point functions T±(q,x) have the same mass-shell 
conditions as in equilibrium, with an off-equilibrium deviation in the spectral densities . 



Furthermore, let us represent the function I^ v (p,x) — / M „ in the form of an expansion 



(7.4) 



where 



u ll {p)u v {p) 
u 2 (p) 



D„ 



are longitudinal and transverse projections in the Lorentz covariant form, respectively, and we take the functions 
in the form of the quasiparticle approximation 



I l p = I l p S(p° - J p ) + lL p S(p° + J p ), P p = I P 5( P ° - w * ) + PL p 6(p° + u p ). 



(7.5) 



In Eq. ( |7.5| ) are certain functions of a wave vector p, and oj 1 ^ = a/'*(p) are frequencies of the gluonic longitudinal 
and transverse eigenwaves in a QGP. 



The equations describing the variation of the spectral densities T q and I l p are obtained from Eqs. (3.1) and ( |6.3| ) 
by the replacements 

T(q, x) -» Mq)T q <5(q° - «-), I^p, x) - Q^(p)[l l p 6(p° - u l p ) + li p 6(p° + u> l p )] (7.6) 
for Eq. (|||), and 

I^(p,x) - Q„ v (p)I p S(p° - oj l p ), T(q,x) -» /i_(q)T q %° - w-) + fr + (q)Tl q % + w~) 

for Eq. ([T^). 5 functions enable us to remove integration over q° and p°. Retaining on the right- hand side the terms 
responsible for the process of nonlinear Landau damping only, we have instead of Eqs. (|6.l|) and (6.3) 



-I ^[<2°+|q|+ Rc <*£_(</)] 



<9T: 



^- = 2 ff 2 C7 F T- /dpl£lmT«>(q,p), 



(7.7) 
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where 



T«> (q, p) = tr{ (ffW) ^(p, -p; q, -q) - T<« "(p; g - p, -g) *S(g - p) TW> "(-p; g, -g + p) 



T W) "(-p; g + P, -g) *S(g + p) T«> "(p; g, -g - p)) h_ (q)}Q^(p) 



(7.8) 



and 



where 



^[p 2 +Re5tf(p)]^ o ^ = -2^4 |dqT q Imr( G )(q,p), 
T( G )(q, p) = tr{(ff (°) ""(p, -p; g, -g) + *r«> "(-p;p - g, g)*S(p - g)*r( G ) "(p; -p + q,-q) 



(7.9) 



T^fa g, -p - g) *S(p + g) "(-p;p + g, -«)) M4)}Q/«<(p) 



The expressions T^(q, p) and T' G )(q, p) contain the factors 

1 1 



(7.10) 



V ■ q + ie ' v ■ g — ie ' u • (g — p) + ie ' 

The imaginary parts of the first two factors should be set equal to zero because they are connected with the linear 
Land au d amping of a plasmino (which is absent in a QGP) and therefore, the imaginary part of the expressions ( |7.8| ) 
and (7.10) will be defined as 



Im • 



v ■ (q~ p) + ie 



-7tS(uj q - ul - v- (q-p)),. 



Because of the absence of th e line ar L andau damping for the plasmino from the de finit ions o f effec tive two-quark-two- 
gluon vertex functions ( |4.7j ) and ( |6.4[ ) and two-quark-one-gluon vertex functions (4.6) and (4.12), the next equalities 
are as follows 



ST (Q) ^(p s _ p . ?) _ q) = ST (G) ^ (pj _ p . Qi * r W) M(p. q _ V: _ q ) = * r ( G ^(p; q - p, -q), 
and as consequence, we have 

TW)(q,p)=r( G )(q,p). 



(7.11) 



(7.12) 



Therefore, in th e sub seque nt discussion the superscripts (Q) and (G) will be dropped. Let us emphasize that the 
equalities (7.11) and (7.12) hold only for on-shell plasma excitations. 



VIII. GAUGE INVARIANCE OF Im T(q, p) 

The issue of the gauge dependence of the plasmino nonlinear Landau damping rate is closely associated with the 
issue of the gauge dependence of ImT(q, p). In the next section the expression ImT(q, p) is shown to define the 
plasmino-hard-particle scattering matrix element squared, i.e., it has direct physical relevance. To establish the gauge 
invariance of the matrix element for plasmino-hard-particle scattering we need to show that the expression ImT(q, p) 
in the covariant gauge equals a similar one in the temporal gauge. 

We prove the gauge invariance for a more general expression, which in a covariant gauge has the form 

*f(pi,p 2 ;gi,g 2 ) = {<5r^(pi,p 2 ;gi,g 2 ) - *r^(pi;gi +p 2 ,g 2 ) *S(qi +p 2 ) *r y (p 2 ;gi,g 2 +pi) 



*r"(p 2 ;gi +pi,g 2 ) *S(qi + pi) *T f2 {p 1 ; gi, q 2 + pa) WfcOM^) 



on— shell 



(8.1) 
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Here, pi + p 2 + qi + q 2 = 0. Notice that the order of the Lorentz indices of t he la st two term s in curly brackets on the 
right-hand side of Eq. (3T) is important. The association of the expression (7J5) with ( |S.l[ ) is given by 



ImT(q,p) 



^u 2 {p)u 2 {pi) 



Imtr|*f(pi,p2;t?i)?2)^-(qi)} 



pi=— P2=p, qi=-q-2=q 



A similar expression (8.1) in the temporal gauge is obtained with the replacements 

"m(p) «/i(p) = P 2 (Pt* ~ u h(p ■ U ))/(P ■ u ), etc. 

The gauge-invariance proof is based on using identities analogous to the effective Ward identities in hot gauge 
theory || M. It can be shown that the following equalities hold: 

5Y^{pi,p 2 ;qi,q2)pi^ = *T u (p 2 ;q 1 ,q 2 +p 1 ) - *T"(p 2 ;qi + pi,q 2 ), 



5T^ u (pi,p 2 ;q u q 2 )p 2v = *T^(p 1 ;q 1 ,q 2 +p 2 ) - *T> 1 (pi;q 1 +p 2 ,q 2 ) 
*Y^p-q 1 ,q 2 )p ll = *S- 1 (q 1 )-*S- 1 (q 1 +p). 



■2) 



Initially we calculate the convolution with the effective two-quark-two-gluon vertex function ST^". Slightly cum- 
bersome but not complicated computations by using the effective Ward identities (3.5) and the mass-shell conditions 
lead to the following expression 



5T'"'(p 1 ,p2; qi,q2,)u M (pi)u v (p2) 



onshcll 



PiP 2 ^T (pi,p 2 ;qi,q 2 ) + E(pi,p 2 ;qi,q2), 



(8.3) 



where 



3(Pi)Pa; 91,92) ^PiPl[*^°(p2;qi +PuQa) ~ *r°(p 2 ; 91, 92 + Pi)j 

+J&?{*r°(p 1 ; 91 + Pa, 92) - *r°(pi; 91, 92 +P2)} - plp^S-'iqi + P2 ) + *S~ x {qi +p x )). 

Furthermore, we calculate the convolution with the terms containing the two-quark-one-gluon effective vertex in 
Eq. (8.1). Here, we derive the expression 

{ *r"(pi;gi +P2,92)*5"(9i + P2)*^ V {P2\ 9i, 92 +pi) + (pi <-> P2, P <-> v) \u^{pi)u v {jp 2 ) 

\. J on — shell 

= PiP2{* r °(Pi;9i +P2,92)*5*(9i +P2)*r°(p 2 ;9i,92 +Pi) + (Pi <->P2)| + 3(pi,Pa; 91,92)- 
Subtracting Eq. ( |3.4| ) from (8.3), we arrive at the desired expression 

*f (pi,p 2l qi,q 2 ) = p\ptiST m {p 1: p 2 ;q 1 ,q 2 ) - *Y Q (j> x ;q x +p 2 ,q 2 )*S{q 1 +p 2 ) *T°(p 2 ; 91, 92 +Pi) 



(8.4) 



T (p 2 ;9i +Pi,92)*S'(9i +Pi)* r °(Pi;9i,92 + p 2 )j. 



(8.5) 



Now we consider the structure of *T in the temporal gauge. For this purpose we replace u M by in Eq. (8.1). The 
convolution with an effective two-quark-two-gluon vertex function leads to 



5T^ (pi , p 2 ; qx , q 2 ) u M (pi )u v {p 2 ) 



onshcll 



PiP2^ r (Pi,P25 9i,92) +3(p!,p 2 ; 91,92), 



(8.6) 



where 



3(Pl,P2!9l,92) = ^^|^Pi*r°(pi;gi +p 2 ,92) - *S 1 {q 1 +p 2 fj + {pi <-> p 2 )j 
Convolution with the terms containing *T M yields 

{*r''(pi;9i +p 2 ,92)*5'(9i +P2)*r ly (p 2 ;9i,92 +Pi) + (pi ^ P2, P ^ ^)}w /1 (pi)5«/(p2) 



shell 



•7) 



2 2 
P1P2 



{*r°(pi;gi +P2,92)*S'(9i +P2)*r°(p 2 ;9i,92 +Pi) + (pi <-> P2)} + S(pi,pa; 91,92)- 



Subtracting Eq. (3.7) from (3.6), we obtain a similar expression (8.5). Thus, we have shown that at least in the class 
of covariant and temporal gauges, the expression ImT(q, p) is gauge invariant. 
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IX. NONLINEAR LANDAU DAMPING RATE OF A PLASMINO 

In our work pl| the transformation of the nonlinear Landau damping rate for a plasmon in a purely gauge sector 
of soft QGP excitations was performed. An expression similar to ImT(q, p) was transformed to module squared of a 
sum of the terms, which are interpreted as the scattering amplitude of specific physical processes. In this section we 
represent a similar transformation for the function ImT(p,q). Since this function is not dependent on the choice of 
a gauge, we choose a temporal gauge for simplicity. 

The right-hand side of Eq. (7.7) contains the contribution of two different processes. The first is associated with 
absorption of a plasmino by QGP particles with frequency uj~ and wave vector q, and consequent radiation of a 

plasmon with frequency wj, and wave vector p. It is defined by the second term on the right-hand side of Eq. ( |7.fi| ). 
The frequencies and wave vectors of an incident plasmino and a recoil plasmon satisfy the conservation law 



v • (q 



0. 



(9.1) 



The second process represents simultaneous radiation (or absorption) of a plasmino and plasmon with frequencies <x> q . 
uj l p and wave vectors q, p satisfying the conservation law 



v- (q + p) =0, 



(9.2) 



and is defined by the third term in Eq. (7.8). The cont ribu tion of the second process to the order of interest is not 
important. The first term on the right-hand side of Eq. ( (7.8|) , associated with an effective two-quark-two-gluon vertex 
function, contains both processes. Fur thermore, we take into account the term s in Eq. (9.1) and we drop the terms 
that contain the S function of Eq. (9.2). Notice that the conservation law (3.1) defines, in a very nontrivial manner, 
accessible kinematic regions of wave vectors q and p in the process of nonlinear Landau damping for a plasmino. In 
Sec. XI we consider only the simplest limiting case of accessible kinematic region of plasmon wave vector p for the 
case of the q = mode. 

Now wc consider the expression with ST^ . With regard to the discussion above and to the definition ( |6.4[ ), the 
contribution of a given term to ImT(q, p) can be represented as 



2lT LOn 



p 2 fdQ (l + vq)(vp) 2 _ , 
pWi^ *K-«P-V'(q-P)) 



(9.3) 



on— she 



Let us consider now a mor e co mplicated second term in Eq. (7^8) associated with the two-quark-one-gluon vertices. 
By the definition (4.5) and (4.6), the following equality is obeyed: 



(hereafter, I = q — p). By using this relation and the decomposition of the effective quark propagator *S(l) onto h±(l), 



where 



*s(0 = Mi)*A+(0 + Mi)*A-(0, 



*A±(0 



Z°T[|l|+<5£±(0]' 



(9.4) 



the contribution to ImT(q, p) of the term associated with the effective two-quark-one-gluon vertex functions can be 
represented as 



P 



P 2 (P°) 



— im(*A + (0tr{/ l _(q)('r i P >+(i)(*rV')}+ *A_(i)tr{fc_(q)(*rV)fe-(i)(*rV)}) 



(9.5) 



on— shell 



Let us compute the traces in the last expression. For this purpose we research the matrix structure of the function 
*r l p l in more detail. From the analysis of the corresponding expressions for HTL functions, derived by Frenkel and 
Taylor (in particular, the expression (3.38) in Ref. |^|), it is easy to see that *T l p l can be represented in the form of 
the expansion 



TV = 7°<5r + (1 • 7 ) T,| + ((n x 1) • 7 ) */j 



(9.6) 
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where n = q x p and the coefficient functions are defined as 

' dfi v ■ p 



Ti, = 



p ' 1 j_ Ar p 1 

~p- + Sr \\ = -p~ 



Uj. 



dfi (v • p) (v • 1) 



Ti = 



l 2 J 47r (vl + ie)(v ■ q) ' 
2 r dn (v • p) (v (n x 1)) 



(9.7) 



l 2 ±_ 1 2 l 2 n 2 7 47r (vl + ie)(vq) ' 



The matrix basis in the expansion (9.6) is convenient in that it is "ortogonal" in trace computing. Substituting 
expression (9.6) into (9.5) we can compute desired traces in Eq. (9.5) in terms of the functions ( |9 . T[ ) . However, this 
direct approach is not quite convenient in view of its nontransparency and the necessity of cumbersome calculations. 
We overdetermine the expansion (9.6), whic h enables us to represent the function T'p* in a more appropriate form for 
trace computing. Instead of expansion ( |9.6| ) we can write (the details of overdetermination are given in Appendix A ) 



TV = - *r+ - Mi) *r- + 2Mq) i 2 |q| *M 



where 



T 4 



-sr oT \n 



|q| lTQ-l 



Ti 



(9.8) 



(9.9) 



The explicit selectivity of the matrix h- (q) in the expansion ( |9.8| ) is connected with the fact that we have restricted 
our consideration to the study of a plasmino branch of fcrmion excitations only. In the c ase of the branch describing 
normal-particle excitations, it is necessary to use the following expansion instead of Eq. (3 

*ry - - m*) % - Mi) *r- - 2Mq) i 2 |q| r±, 

where now 



(9.10) 



t± = - sr T |1| T 



1 



|q| l±q l 



T, 



(9.11) 



Substituting expression (9.8) into ( |9.5[ ), and using the identities for h±(l) matrices, 

h±(i)h ± (l) = o, Mi)Mi)Mi) = h ± (l), 

after computing the trivial traces, we obtain 



(1 - q- 1) Im (*A+(0(T+) 2 ) + (1 + q • 1) Im (*A_(/)(T_ 



on— shell 



p 2 (p ) 2 

Furthermore, we use the relation 

Im(*A±(/)(T ± ) 2 ) = -Im(*A ± 1 (/))|*A±(/)T±| 2 + 2Im(T±) Re(*A±(/) T±). 
Taking into consideration the last relation and the equalities 

Im *A^(0 = - ™lj ^ (1 rp v • i ) S(v ■ I), 



(9.12) 



ImT± 



lTv-1- 

47T v ■ q 



v • (n x 1) 
|q|l 2 (lTq-i; 



6(vl), 



subtracting Eq. (9.12) from (9.3), we define the desired expression for ImT 



ImT(q,p) 



7TW n 



p 2 (p ) 2 



d^L^p) 2 

4tt \ (v q) 2 V 



v- q , 
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(l-q-l)(l-vl)|*A + (0*r + | 2 + 2 



v P 

v ■ q 



2 + (v;q,p)Rc(*A+(0T + ) 



(9.13) 



(l + q-l)(l+vl)|*A_(0*r_ 



v P 

v ■ q 



e _(v;4,WRe(*A_(0*r_)U(vO 



- on — shell 



Here. 



g±(v;q,p)E(l T q.l)(l T v.i; 



v ■ (n x 1) 

|q|i 2 ' 



(9.14) 



The expression ( |9.1£ ) can be led to a more descriptive form. For this purpose we add to the expression in the curly 
brackets on the right-hand side of Eq. ( |9.13[ ) the term 



v • (n x 1) 
|q|l 2 



(rA+(Z)T + | a + |*A_(0T_| a ), 



equals zero over solid integration. In addition, we rewrite the first term in curly brackets by using the relation 

2(1 +vq) = p + (v;q,p) + £>-(v;q,p), 



which is a consequence of the definitions (9.14). Taking into account the discussion above, we rewrite expression 
( |9.13 ) in the following form: 



I \2 



ImT(q,p) 

where 



:K) 



K) 2 p : 



p2 /> 

J ~^ s ^ _CJ p ~ v (<J - p)){?+( v i4p)^(qiP) + e-(v;q,p)«v (i.p)}i ( 9 - 15 ) 



Wv ± (q,p) = |X ± (q,p)|^ = 



v P 

V ■ q 



+ *A ± (z)*r ± 



(9.16) 



on — shell 



The denominator v ■ q in Eq. ( 9.16 ) is eikonal, which was expected in an approximation to the small-angle scattering 
of a high-energy particle. It defines the effective (here, effectiveness is not in terms of HTL resummation) propagator 
in the Compton type scattering process depicted in Fig. 4. The factor v • p is connected with the effective vertex 
with one soft external leg and two hard ones. Accordingly, the second term in the amplitudes A4±(q, p) represents 
the second type of scattering process depicted in Fig. 5. Here, the quark propagator transfers the soft momentum 
and is defined by the scalar functions *A + or *A_, and a vertex with soft external momenta is defined by the scalar 
functions *T+ or *T_ for corresponding propagators. 

Furthermore, we are relating the Wigner functions Tq and I 1 to the plasmino and plasmon number densities, 
setting accordingly 



Here. 



z: 1 (q) = i 



q =2ZI 1 (q)T q , Nl = -2j p Z^(p)I l p . 

fdReSU l (p) 



(9.17) 



,o 



ZP(P) 



1 



V d( P °) 



are the residues of the effective quark and gluon propagators at the appropriate poles, respectively. The factor 2 in front 
of Z I 1 (q) takes into account the presence in QGP's of antiplasminos. By using the definition ( 9.17 ), Eqs. (7.7) and 
(7.9) can be rewritten in the usual form of kinetic equations (containing on the left-hand side drift terms, and on the 
right-hand side the terms responsible for collisions) , where the role distribution functions of quasiparticles-plasminos 
and plasmons-fulfill the functions (9.17) 



at 



V; 



5x 



= - 7 (q)n = - g 2 C F n / dp Q(q, p)JV' 



(9.18) 



<9x 



-7'(PK 



-g 2 T F N l p Jdc i Q(q 1 p)n^ 1 



(9.19) 
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where V q = 9w q / 9q and V p 
<2(q, p) is denned by 



dui 1 I dp are group velocities of plasminos and plasmons respectively, and the kernel 



Q(q,P) = ^ r Z_(q)Z i (p)ImT(q,p). 
2w p 



(9.20) 



The functions 7~(q) and 7 Z (p) on the right-hand side of Eqs. ( 9.1 8| ) and (9. IE) represent the nonlinear Landau damping 
rates for plasminos with mome ntum q and for plasmons with momentum p, respecti vely. 

The structure of kernel ( 9.20| ) is rather unexpected. As we see from the expression (|9.15|) , this kernel is not reduced 
to the squared modulus of one scalar function, as this occurs in a purely gauge case |21||. Here, Q(q, p) is defined by 
a sum of the squared moduli of two independent scalar functions 5 : amplitudes .M + (q, p) and A4-(q, p). This is a 
point that we find difficult to interpret and therefore additional analysis of this problem is required. The only remark 
that may be made is that the interference between the scattering processes depicted in Fig. 5 proceeding through 
intermediate quark virtual oscillations with propagators *A + (Z) and *A_(Z) accordingly, vanishes by the relation 



Mi) 



■<4>-v(q-p)) 



Mi) = o. 



This relation is the analog of the relation in a purely gauge case [ElJ 



(Pa) 4 



J ^_ V V 5{J p - u l pi - v • (p - pi)) 



(n x p 2 y =0, n = p X pi, p 2 



P - Pi, 



that is responsible for the absence of interference between the scattering of plasmon by a QGP thermal particle 
through the longitudinal and transverse virtual gluon oscillations. 

X. DECOMPOSITION OF KERNEL Q(q,p) INTO POSITIVE AND NEGATIVE PARTS 



In this section we consider the problem of the direction of the effective pumping over of plasma excitation energy 
in the process of the nonlinear interaction of plasminos and plasmons. For this purpose, first, we study in more detail 
the structure of the g±(v; q, p) function. 



By using the expansion (A4), we represent the functions g± (9.14) in the form 



e±(v;q,p) = l + v-q=F (q-l + v-1). 



(10.1) 



Let us introduce the coordinate system in which axis 0Z is aligned with vector 1; then the coordinates of vectors q and 
v are equal to q = (1, a, f3) and v = (1, 0, <p), respectively. By <I> we denote the angle between v and q: v • q = cos $. 
The angle $ can be expressed as 



cos $ = sin 9 sin a cos(<y9 — j3) + cos 9 cos a. 



In the fixed coordinate system expression (10.1) reads 

g±(a, (3;9,ip) = 1 + sin 9 sin a cos(<p — (3) + cos 9 cos a =F (cos a + cos ( 



(10.2) 



(10.3) 



By using trigonometric formulas, it is easily to show that the last expression can be represented in the homogeneous 
quadratic form with respect to the variables cos |(0 — a) and cos |(0 + a) 



Q±(a,f3;6,(p) = yU±y 



(10.4) 



where 



y(M) 



cosi(#-a)\ . . , R /2cos 2 i((^-/3) T l 

C o4\e + a))' A±=A± ^=\ T l 2sin 2 I(^-/3) 



5 As will be shown in the next section, the coefficient functions g± in Eq. (9.15) are not, in general, case definite; therefore in principle it 
is impossible to express the expression Im T as the squared modulus of one scalar function-the total scattering amplitude. 
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The eigenvalues of matrices A± are equal to 

A (+) = 1 + yj 1 + cos% - J) > 0, A ( ~ } = 1 - a/T+ cos 2 (v? - /3) < 0. 



(10.5) 



From the theory of matrices p8fl it is known that the real symmetric qua dratic form is indefinite if and only if the 
eigenvalues A of matrix A have the different signs. Thus, by virtue of Eq. (10.5) we have proved that the coefficient 
functi ons g ± are indefinite. By the linear transformation of the variables cos ^ ifi ~ Qf) and cos ^(8 + a), the quadratic 
form (10.4) can be expressed as the canonical one, where the positive and negative parts of g± are explicitly displayed 



Q± 



(a,/3;0,^)=A( + )(xL +) ) 2 + A(-)( X i- ) ) 2 . 



(10.6) 



Here, x± an d x± are certain linear functions of cos ^(8 — a) and cos ^(8 + a). 

The explicit expressions for x±^ an d X± are defined with the help of the transformations 



xL +) 



= u ± 



cos \ (8 — a) 
cos 



(10.7) 



> yi> ~ \~ ck) 

where, in turn, the transformation matrices U± arc defined by the solutions of the corresponding matrix equations 

U±A± = diag(A (+ >,A ( -))W±. 
The solutions of the last equations have the following structure 



W± = 



cos i? ± sin d 
=F sin ■& cos d 



teai-d = — 



1 



(<P-P) + Vl + cos 2 ((^-/3)' 



cos 



By using a given expression and Eq. (10.7), we derive the desired expressions for x±^ an d X± The functions Xdt"^ 
and x± "* can be easily rewritten in the terms of initial vectors 1, q and v. For example, from Eq. ( 10. 2j ) it follows that 



cos(ip — 0) 



n • (v x 1) 
|n||v x 1| 



etc. 



By substituting the expansion ( 10.6 ) into ( 9.15 ), we can represent the kernel Q(q, p) in the form of a sum of positive 
and negative parts 



Q(q,p) = Q (+) (q,p) - Q H (q,p), S (±) (q,p) > o, 



where 



S (±) (q,p) = ±™ 2 



K) 2 -p 2 



2K) 3 P 2 



Z-(q)Z,(p) 



(10.8) 



(10.9) 



47T 



v ■ (q - P) ) { A<±) [(xi ±} ) 2 ™+ (q, P) + (X W ) 2 «V (q, P)] } ■ 



We use the decomposition (10.8) in general analysis of the problem of a direction of the pumping over of excitation 
energy in the study of the process of nonlinear interaction between plasminos and plasmons. 

To consider the model problem of interaction of two infinitely narrow packets with typical wave vectors qo and po, 



let us introduce the number densities n q and as follows: 



n-(t) = Ti- (t) *(q - qo), NL(t) - N l (t) 6(p - p ). 



We h ave restricted ourselves to the spatially homogeneous case. Substituting the last expressions into Eqs. (9.18) and 
we obtain the coupled nonlinear equations: 



dn 
~&t 



= -g z C F Qn-N\ n-(t ) 



"o ■ 



(10.10) 
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8N l 

= +g 2 T F Qn-N\ N l (t ) = N l . 

Here, Q = Q(qo, Po)- The system of the equations ( |10.10D possesses the integral of motion 

C = T F n-(t) + C F N l (t) = T F n +C F N l . 



The general solution of this system, with regard to decomposition of the kernel ( 10.8 ), has the form 

e -ff 2 CQ (+) (t-t ) 

U (t) = n ° ° C F N l e-9 2 CQt-)(t-t ) + T F n c-9 2 CQ^(t~t )> 

e -9 2 CQ<-»(t-t ) 

W C F N l e-9 2 CQi-)(t-t Q ) + Tp n - e ~9 2 c Q (+) (t-to) • 

Let us analyze the behavior of these solutions in the limit for t — > oo. 

1. Let us assume that the values of the wave vectors qo and po are such that the following inequality 

Q (+) (q ,Po)> G (_) (qo,Po) (10-11) 

is true. Then in the limit for t — ► oo we have 

n-(t)-fO, N\t)^ N l + {T F /C F )riQ. (10.12) 

Thus we see that as a result of the nonlinear interaction of two infinitely narrow packets with fermion and boson 
quantum numbers, the effective pumping over of energy from the first packet to the second one takes place. 

2. If the inverse inequality holds, 

Q (+) (qo,Po) < Q ( " 3 (qo,Po), (10.13) 

then in the limit for t — > oo we have 

n-{t)^n- + {C F /T F )Nl N l (t)^0. 
Here, we have an inverse case-the damping of bosonic excitations and the growth of fermionic excitations. 

3. Finally, in the limiting case Q^ + '(qo,Po) = '(<10)Po)j the interaction between wave packets is absent. 



The inequalities (10.11) and (10.13) define the kinematic relations between wave vectors of excitations with different 
statistics, such that one or another process of pumping over of energy occurs. However, the general analysis of these 
inequalities is a complicated problem by virtue of the complexity of the expressions for positive and negative parts of 
the kernel Q. In the next section we shall study in detail only the simplest limiting case of the interaction of "standing" 
plasmino (q = 0) with a plasmon. Here, the direction of the effective pumping over energy is easily defined. 

XI. THE DAMPING RATE OF A PLASMINO AT REST 



The solution of the kinetic equation (9.18), defining a change of the plasmino number density, can be formally 



represented in the form (for a spatially homogeneous case) 

t 

*4 -) (t) = < _) (*o) ex P {- g 2 C F J dt'J dp Q(q, p)<(*')}- 

to 

Let us approximate the plasmon number density by its equilibrium value or Planck distribution 

<-< (|p|) =(2^e^TT 
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and define the quasiparticle damping rate of the standing plasmino by means of the relation 

7"(0) = \g 2 C F lim I dpQ(q, P X(|p|). 



(11.1) 



Now we represent a complete calculation of 7~(0). We start wit h a r eprese ntatio n of a kernel <2(q, p) in the form 
(5). 20), where the function ImT(q,p) is determined through Eqs. ( |9.15 ) and (9.16). Here, it is convenient for us to 
introduce a new coordinate system in which axis OZ is aligned with the vector p; then the coordinates of vectors q 
and v are equal to q = (1, a', (3') and v = (1, 6', (p 1 ), respecti vely. By $' we denote the angle between v, and q. The 
angle $' is connected with a' , /?', 6' and p' via the relation (10.2) with a corresponding change a — > a', etc. 

In the limit for |q| — > 0, the kernel Q(q, p) is reduced to 



Q(0, p) = ™l K of P ;, ^ HP) I ^ S(cose' - p l p ){ g+ (v- q, p)|M + (0, p)| 2 + g_(v; q, p)|M_(0, p)| 2 }, (11.2) 



K)3|p|3 



where p p = (u> p — wo)/|p|; dfl = sin 9 ' AO' dip' . The expressions for £?±(v;q, p) are defined by Eq. ( 10.1 ). In a new 
coordinate system, instead of Eq. ( |10.3| ) , we have the function representation g± in terms of angles 

Q±(a',(3';6', tp') = l + sin 6»' sin a' cos (y/ - (3') + cos 6' cos a' ± (cos a' + cos0')- (H-3) 



The limits of the scattering amplitudes Ai±, by its definitions (9.16) equals 



M±(0,p) = i^cos6»'+ lim (*A±(n*r±). 

UJ |qH0 



(11.4) 



Thus, the problem of computing of 7~(0) is reduced to the calculation of limits of scalar functions: the effective 
propagators *A± and the effective vertices *T±. 

By using the definitions of the *T± effective scalar vertex functions Q.S) and ( |9.7| ), after slightly cumbersome 
computations, we define 



lim T± 

|q|-»0 



IPI J 



■ p» p lim *A ± X (0 - [ co T | P | T — P p (l ± Pp) ] • 

LU y |q|—0 W ^ P 



(11.5) 



In the last eq uality we us e the *A ± 1 (g) definition (3.4) as the F(go/|q|) function (7.2) (see Braaten and Pisarski (g[). 
Inserting Eq. ( 11.5 ) into ( 11.4 ), we reduce the scattering amplitudes to 



M±(0, p) = M cos0' - — p p - [w =F |p| T — P P (1 ± P P ) ] , lim *A±(Z). 



By the S function in Eq. (11.2), all terms in A4±(0, p) not containing the propagators *A±, are relatively reduced in 
the limit of the q = mode. 



The remaining terms, after substitution into Eq. (11.2) and integration over the solid angle with regard to Eq. (11.3) 
yield 



Q(o,p) = 



K) 



i ^2 p2 



16(w p ) 3 |pl 3 



Z,(p)0(l - 1^ |) (l + cosa')(l + P P ) Wo-|p|-— Pi(l + P p ) lim*A + (0 



OJQ 



Iqho 



+ (l-cosa')(l-P P )[^o + |p| + — P l p (l-P p )] 2 lim "A_(0 
L - 1 |q|-*o 



(11.6) 



We note that this expression i s not dependent on the angle f3' . This enables us to represent the integration measure 
on the right-hand side of Eq. (11.1) in the form 



oc — 1 

ydp = 2nJ p 2 d|p| J d(cosa'). 



Substituting Eq. (11.6) into (11.1) and performing an angular integration over a' , we finally obtain 7 (0) 

00 

7"(0) = g 2 C F J d|p| 0(1 \p l p \) e-(|p|)^(|p|), 



(11.7) 
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where the kernel Q (|p|) has the form 



Q-(|p|) = ^ 2 



H) 2 -p 2 

4K)3|p| 



p 2 j 



z,(p)ui+p') wo-ipI-— Pp (i+p; 



w 



lim *A+(l) 

|q|-»0 



(11. 



P P ) 



w 



-Pp(l-Pp) 



lim *A_(/) 

|q|->0 



By virtue of the 9 function in integrand ( |ll.7| ), the kernel Q (|p|) is positive; therefore 

7~(0)>0, 

i.e., the standing plasmino is damped. Thus, in the process of nonlinear Landau damping, the pumping over of the 
excitat ion energy of standing plasminos into plasmon branch of plasma excitations occurs, and therefore the first case 
(|10.12|) is true. 



The function p , entering the argument of 9 function and in a kernel ( 11.8 ), decreases with momentum |p| from 



+00 reaching a minimum, and then monotonically increases at large |p|, asymptotically tending to +1 from below. 
The equation p = 1 defines a lower limit of integration over the momentum of the recoil plasmon. The numerical 
solution of a given equation yields 



0.447 u p i for N f = 2 



/ 



0.495 w p i for N f = 3 



(hereafter, the numerical estimates are presented for the SU(3) color group). The function p l p reaches an absolute 
minimum at 



1.236 u) 



'pi, P P ** 



0.664 for N 



f 



1.305 Upi, p p ,, ~ 0.692 for Nf = 3. 



Notice that the value of the function p l p at a minimum point coincides with the value of the plasmon group velocity 
do>p/d|p|. Thus, if the incident plasmino is initially at rest, both the recoil plasmon and the effective quark that 



propagates below the light cone carry nonzero energies and momenta of order gT. 



of standing 



At the end of this section we compare derived expressions for the damping rate ( |11.7|}_ and ( 11. 
plasmino with the similar expression obtained in the framework of the HTL approximation ]|, ^] . For this purpose in 

Wq|p|- The Planck distribution JV' (|p|) in the integrand in 



Eqs. (11.7) and ( 



11. 



we rescale uj, 



■* uj uj p and |p| 



Eq. (11.7) should be set equal T/((27r) a o;o S p ), since the energy is soft. Further rewriting the kernel ( 11.8| ) in terms 
of the functions accepted in the paper || and defining the damping constant as 



j-(0) = a(N c ,N f ) 



9 2 TC F 
4tt ' 



we derive the expression for the coefficient a(N c ,Nf) required for comparison 



a(N c ,Nf) 



IpNpI 



|p|)0(|p|-|5' -1|) 



:{[ |p| + - 2 ] 2 /?+(l - 2' , |P|) + [ |p| - 2 p + 2] 2 /3_(l - w' |p|)}, 



where 



Ri(po, IpI 



Po Pp - P 2 
p 2 3r - pi + p 2 



For the definitions of more cumbersome functions /§±(po, |p|) see Ref. The coefficient a(N c , Nf) exactly coincides 
with the corresponding part in a similar coefficient a(N c , Nf) in the expression of the damping rate of the standing 
plasmino, responsible for the scattering process of plasminos by hard-particle QGP's varying with the change of 
excitation statistic, derived by Kobes, Kunstatter, and Mak ( |9]|, E q. (5.11)). We have shown in doing so, by really 
remaining in the context of the Blaizot-Iancu equations (2.1)-(EJj), we are able to compute not only a spectrum of 
fermionic excitations in QGP's, but also their (gauge-invariant) damping rate at the leading order in the coupling 
constant, which corresponds to the damping rate from the resummed perturbation theory [& H|. 
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XII. LIGHT-CONE SINGULARITY: IMPROVED BLAIZOT-IANCU EQUATIONS 

In this section we discuss the behavior of the nonlinear Landau damping rate for normal-particle excitations, when 
the quark spectrum w+ approaches the light cone. Here, the initial expression is 

7 +(q) = i 5 2 C F ydpQ(q,p)A^(p), Q(q,p).i r Z + (q)Z I (p)ImT(q,p) I (12.1) 



where 



{uL ) 2 - p 2 fdfl 

ImT(q,p) = ttu ^ ^ 2p2 J — - u> p - v ■ (q - p)){g+(v; q, p)w+(q, p) + g_(v;q,p)w v (q,p)}, (12.2) 



(t^) 2 p 2 J 47T 

and 



e ± (v;q,p)E(l±q.i)(l T vi)- V ' ( " * 1} , (12.3) 

|q| V 



wt (q>p) 



l-M±(q,p)| 2 



v P 



+ *A±(Z)*/± 



(12.4) 



shell 



The scalar vertex functions *T± are defined by Eq. ( 9.11 ). On the right-hand side of Eq. ( 12.1 ) we have taken into 
account for simplicity only the contribution from the nonlinear interaction with longitudinal bosonic excitations, but 
all subsequent derivations easily extend to the case of the nonlinear interaction with transverse excitations. 

Near the light cone, the HTL- vertex pieces Sr , SFp and ST±, the explicit forms of which are given by Eq. (p. 7]), 
are the origin of the strongest singularity of order 1/e, where e 2 = ((ujt) 2 — q 2 )/q 2 - The terms linear in <5T in the 



probabilities (12.4) can lead to the logarithm of e only |15[| (the first term in the amplitudes A4±, defining Compton- 
like scattering processes, also results only in logarithmic divergence). Hence we restrict our consideration to the terms 
quadratic in Sr. 

Now let us single out in the HTL- vertex pieces the terms generating near- light-cone 1/e-contributions to 7 + (q). 
As an example we consider the scalar function Sr . By using the explicit expression for HTL amplitudes derived by 
Frenkel and Taylor M , we find 



v p 



dn 



47r (v ■ I + ie)(v ■ q) 



2 f dil V v ■ p 

47r v ■ I v ■ q 



~ — [ Z°(n X q) • p - q°(n x 1) • p ] M(Z, q) = - uj 2 oP °M(1, q). 

In the above expression, V indicates the principal-value prescription; here and in the following ~ indicates that we 
have dropped less-divergent terms. The Lorentz-invariant function WL(l,q) is defined by 



I J l-q+^-AQ^ Y A(l,q)^Pql-(l.qr<0 



M(l,q) = { 



2v/-A(Z,g) \l-q- y/-A(l,q) 
1 



■. arctan 



l-q 



(12.5) 



A(l,q)>0. 



Similarly, the terms containing the function M(/, q) in the HTL- vertex pieces 5T|| and 5T± are singled out. Substituting 
next thus derived expressions into Eq. (9.11), we obtain 



1± ~ u 



( q »-n±( q »+n 



o , l(h <Z°|q|Ti°|l| 



|q||l|±q-l 



M(i,g) = wg/ ± (|l|,|q|,l.q)M(i,g). 



Furthermore, substituting the last expression into the probabilities (12.4) and keeping only the terms quadratic 
in M(l,q), we have 



w$(q,p) ~ ^ |*A ± (Z)| 2 / 2 |M(Z,g)| 2 , /± = /±(|1|, |ca|, 1 - q). 
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The depende nce o n the unit vector v in the integrand of the expression ( 12. 2j ) is defin ed only by the 5 function and 
g± f uncti on ( 12. 3| ), that enables easily to perform the angular integration in Eq. ( 12. 2|) and thus to derive instead of 
Eq. fll2~2]) 



ImT(q,p) 



P 2 1 



° K) 2 P 2 2|1| 



9(-l 2 ) 



(12.6) 



: {(1 - jjj)(l + q • i)/ 2 |*A + (0| 2 + (1 + - q ■ l)f_ |*A_(0| 2 } |M(I, q)f 



The function ImT(q, p) is different from zero for I 2 = (q — p) 2 <0. It is not diffic ult t o show that the condition I 2 <0 
on the mass shell of plasma excitations leads to A(/, q) <0, and therefore in Eq. (12.5) it is necessary to take the first 
expression. 

With q 2 — > 0, the square of M(l, q) becomes 1/e times a representation of the 8 function 



\M(l,q)\ 




S(q-l). 



(12.7) 



Substituting further the expression ( 12. 6| ) into Eq. ( 12.1 ) and taking into account Eq. ( 12.7 ) and Z + (q) q — ^ 1, we 



find the most singular contribution to the nonlinear Landau damping rate of the normal quark excitations near the 
light cone: 



- \^C W m / d |p| *(P)AL(P) 



e 2|q| 



■in 



K) 3 |i|v=P 



(12. 



: {(1 - £)(1 + q • l)fl |*A + (0| 2 + (1 + jMl q ■ l)/ 2 |*A_(0| 2 } 9(-l 2 ) S(q ■ I). 



Here, the solid integral is over the directions of the unit vector p. The S function in integrand (12.8) enables us 
in principle to perform the angular integration. Because of the fact that the functions Z/(p) and A^ q (p) vanish 
exponentially at large |p|, the above exp ression is ultraviolet convergent. Notice also that by virtue of the 9 function 
in the integrand, the damping rate (12.5) is positive. Thus we have shown that really the nonlinear Landau damping 
rate of the normal quark mode near the light cone diverges as \j\fq*', and thus signals the need to have further 
improvement of the Blaizot-Iancu equations, beyond the HTL approximation. 

As mentioned in the Introduction, the light-cone singularities are associated with the massless basic constituents 
of a hot plasma, in our case, with massless hard quarks and hard transverse gluons. The inclusion of the asymptotic 
thermal masses for the basic constituents removes this type of singularity. The effective way of entering thermal masses 
into the resummed perturbation theory without spoiling gauge invariance was suggested by Flechsig and Rebhan jl7) . 
Below we use this approach, reformulating it into a "kinetic" language. 

T he lig ht-cone singularity is generated by the HTL piece of the effective two- quark-one- gluon vertex functi on (|4.(i| ) 
(or (|DJ)). In turn, the expression for <5 F(jp; q±, q2) is defined with the help of the Blaizot-Iancu equation ( |2.8| ) for 
the function dn± or by application of Eq. ( [2.9|) for function I)L. To derive the expression for S T(p; q i , (?2 ) , which is free 
from the light-cone singularity, we use improved Blaizot-Iancu equations instead of Eqs. ( |2.8D and ( |2.9| ) (the details of 
deriving these equations are given in Appendix B) 



[v ■ D x ,Snt(k, X)]=±£- t a (4>(X)t a £±(k, X) - A^k, X)t a i,(X)), 



(2.8') 



(v ■ D x )^(k,X) ± #±( k ,X) = - igC F [N(e k ) + n(e k )]^(X), (2.9') 

where J ^ ± (k, X) = (jX ± (k, X)Yj , Am 2 ^ = m 2 q oo — m 2 g oo , and mn x and m g:00 are asymptotic thermal masses for 
hard quarks and gluons, respectively. Instead of the expression (fejs) for the induced source now we have 



V(X) 



lI(^^J + ^ + °Jwf^ k ™ *)• (2.30 
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The expression for the color current induced by the soft-quark field (2.5), is not changed. 

Employing the i mp roved Blaizot-Iancu equations (2.8') and (2.9'), and expressions for induced source (2.3') and 
the color current (2.5), it is not difficult to obtain the improved soft-quark self-energy and the effective two-quark- 
one-gluon vertex function used above. These expres sion s coincide with similar ones derived in the context of thermal 
field theory p7|. Thus, for example, instead of Eq. (4.6) now we have 




qi 



2Ta 



+ te)(v ■ q 2 



2Ta 



(• 



2Ta 



(• 



le \v ■ q 2 



2Ta 




(12.9) 



Here. 



2 f ae a 
(0(a)) a = — da^-jO(a), a = e k /T. 



The use of expression (12.9) gi ves a finite nonlinear Landau damping rate 7 + (q) near the light cone. The additional 
ter ms in denominators of Eq. ( 12.9 ) lead to the replacement of the divergence factor 1/e on the right-hand side of 
Eq-dU): 



T\q\ 

Ami 



1/2 



where (a 1 / 2 ) Q = (3/47r 3 / 2 )(l — 2~ 2 / 3 ) £(5/2) and £ is a Riemann zeta function. 

At the end of this section notice that we run into similar divergence in research of a behavior near the light cone 
of the nonlinear Landau damping rate for the transverse bosonic excitations in purely gauge sector |^lf . The light 
cone singularity here is generated by the HTL piece of the effective three-gluon vertex function. It is clear that in 
this case for elimination of this singularity it is necessary to use the improved Blaizot-Iancu equations for functions 



Sn± and 5N A instead of Eqs. (2.6) and (2.7). Here, we have not presented their explicit forms, noting only that their 



derivation is rather different from deriving improved equations (2.8') and (2.9') 



XIII. CONCLUSION 



In Sec. X the nonlinear interaction of plasminos with plasmons is shown to lead to the effective pumping over of 
plasma excitation energy from the fermionic branch of plasma oscillations to the bosonic branch and vice versa, and 
the kinematic relations between wave vectors of excitations are defined, such that one or another process of pumping 
over of energy occurs. However, it is clear that conclusions made in Sec. X are somewhat restricted since they were 
made with allowance only for one process-the process of nonlinear Landau damping. For precise study of nonlinear 
dynamics of soft-fermion and soft-boson excitations in hot QCD plasma it is necessary also to take into account 
in the plasmino and plasmon kinetic equations the remaining terms on the right-hand sides of generalized kinetic 
equations (6.1) and ( |6.3| ), responsible for the decay processes etc, and consider the existence of s oft no rma l quark 
excitations and soft transverse gluon ones (the kinetic equations for them are defined from Eqs. (6.1) and (|6.3|) by the 
replacements similar to Eq. ( |7.6[ )). In addition to the right-hand side of the generalized kinetic equation ( |6.3D should 
be supplemented with the terms from purely gauge sector |2l[ . 

Thus obtained self-consistent nonlinear system of four kinetic equations will contain the maximum comprehensive 
information on soft excitations dynamics in hot QCD plasma, which may be obtained in the context of the initial 
Blaizot-Iancu equations (2.1 )— (pT3) (or its improved variant) in the first nonvanishing approximation 6 in the expansion 



The following term in the expansion over this small parameter leads to the processes of Boltzmann type: 

q + q(q) ^ q + q(q), q + q^g + g, g + g^g + g, etc., 

going without exchange of energy between hard thermal particles and soft plasma waves. The example of the construction of the 
Boltzmann equation, describing the elastic scattering process of colorless plasmons between themselves can be found in Ref. B9|. 
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over the small parameter gRA^, where R ~ (gT) _1 and |^4 M | <§C T. Certain nonlinear processes will be kinematically 
forbidden and therefore relevant contributions on the right-hand side of kinetic equations will drop out. However, 
even with allowance for the last case, this system of equations is especially complicated for any analytical research 
(similar to studies carried out in Sec. X of the present work and in Sec. 10 in Ref. |2l|), and therefore here invoking 
numerical methods is required. 
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Appendix A 



By using the relations 7 = /i_(l) + h + (V) and 1 • 7 = h-(V) — h + (l), we rewrite the expansion (9.6) in the form 

*ry = Mi) (flo + |i| t,|) + Mi) (sr - |i| T,,) + ((n x i) - 7 ) *r ± . (Ai) 

Let us add to the expression in the first parentheses on the right-hand side of Eq. (Al) and then subtract from 
it the function n 2 /(|q|(l — q • 1)) and in the second parentheses add and subtract the function n 2 /(|q|(l + q • 1)) 
correspondingly. Let us introduce the following functions: 



T ± = -<y + |i|T| 



Instead of Eq. (Al), now we have 

TV = - Mi) * r + - Mi) * r - 



q| Vi-q-i 



|q| lTQ-l 



Mi) 



T, 



1 + ql 



h+(l) +((nxl)- 7 ) 



(A2) 



Furthermore, we transform the expression in the parentheses on the right-hand side of Eq. (A2). By using the 
definitions of matrices h±(l), this expression can be rewritten in the form 



7 



1-ql 1+ql 



(1-7) 



1-q-I 1 + ql 



q2l2 (7° + (i-7)(q-i)) = 2^Mq) 



ri- 



li- 



ng 



((nxl). 7 ). 



(A3) 



In deriving the last equality we have used the expansion of vector q into two mutually ortogonal vectors, 1 and n x 1: 

q = V 1 + T^ (nxl) - (A4) 

Substituting Eq. (A3) into (A2) and reducing the terms with vector product, finally we derive the expression (|9.8j ) 
for convolution *F*p l , instead of Eq. (9.6). 



Appendix B 



Here we show how one can obtain the improved Blaizot-Iancu equations, used in Sec. XII for deriving the finite 
nonlinear Landau damping rate for normal-particle excitations near the light cone. As the initial equations for 
obtaining an improved equation for a one-body density matrix mixing fermion and boson degrees of freedom, we 
take Eqs. (3.33) and (3.34) in the original paper of Blaizot and Iancu |22] (see the accepted notations and definitions 
therein) 



[d 2 s +d s -d x + 2igA(X) ■ d s ] KZ<(s, X) = ig (^A) 7 Tf(I), 
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KV(d s 2 - d s ■ d x ) ~ gf abc T^ pX AP(X) c\ A ] K»<{s,X) = - gS<(s)^(X). (B2) 

Furthermore, we expand K< = + ..., with if W ~ gK {0 \ etc. The dominant terms in Eqs. (Bl) and (B2) 

lead to the consistency condition 

d 2 s K {a \s.X) = 0. (B3) 

Equation (B3) can be considered as the Klein-Gordon equation for hard, massless, free particles. However, due to 
interaction with the hot medium, the basic constituents of the plasma are known to acquire the dynamical (asymptotic) 
temperature-induced masses. To account for this fact we modify the condition (B3), replacing it in the case of Eq. (Bl) 
by 

d 2 s K(°\s,X)^-ml 00 K^(s,X), (B4) 

where m 2 ^ = g 2 T 2 Cf/4 is the asymptotic thermal mass for hard quark, and in the case of Eq. (B2) by 

d 2 s K^{ S ,X)^-ml 00 K^{s,X), (B5) 

where m 2 ^ = g 2 T 2 (Nf + 2N C )/12 is the asymptotic thermal mass for the hard transverse gluon. It is necessary to 
keep the terms d 2 inside the brackets on the left-hand sides of Eqs. (Bl) and (B2) setting them equals the right-hand 
side of Eqs. (B4) and (B5), accordingly. Further, subtracting Eq. (B2) from (Bl) and performing transformations 
similar to those in Ref. p2[ , we come to improved equation for a density ^Z(fc, X) = tf^^K^k, X): 

(v ■ D X )$%X) - -Am^ ${k,X) = -igC FPQ (k) [N(k ) + n(k )} )6 i/j(X). 
Here, Am^, = m 2 x — m 2 00 . Substituting into the last equation the expansion 

#(fc,JC) = 27n5(fc 2 ){(9(fc )A /+ (k,X) + 0(-fc o )4T(-k,Jr)}, (B6) 

replacing k — > — k, where it is necessary, and dividing by we obtain the improved Eq. (2.9'). By virtue of the 
fact that Am 2 ^ ^ 0, the discrepancy ^ + (k, X) ^/L~(k, X) will hold. The last case is an actual reflection at the 
kinetic description level of the discrepancy of the two terms in the integrands of the improved hard thermal loops 
with external fermions when exchanging k — > — k, found by Flechsig and Rebhan |l7j . 

Now we turn to deriving the improved kinetic equations for <5n*. The equations for functions <5n* are defined from 
the equation for the current density J*(fc, X) ]22Tj : 

[v ■ D x ,J*(k,X)} = tgk^t a {4,(X)t a ^(k,X)- nk,X)t a ^j(X)} (B7) 

by substituting to it the expansion J^(k.X) = k^4:Trd(k 2 ){8(ko)Sn^_(k, X) + 9{— fc )5n*(— k, X)}. The relevant 
modification of the equation for Sn± here is achieved by the simple requirement of using the expansion (B6) on the 
right-hand side of Eq. (B7), where the functions ^ obey the improved equations (2.8'). 



[1] V.V. Klimov, Yad. Fiz. 33, 1734 (1981) [Sov. J. Nucl. Phys. 33, 934 (1981)]; Zh. Eksp. Teor. Fiz. 82, 336 (1982) [Sov. 

Phys. JETP 55, 199 (1982)]; H.A. Weldon, Phys. Rev. D 26, 2789 (1982). 
[2] R.D. Pisarski, Nucl. Phys. A498, 423c (1989). 

[3] E. Braaten, Astrophys. J. 392, 70 (1992); G. Baym, J.-P. Blaizot, and B. Svetitsky, Phys. Rev. D 46, 4043 (1992). 
[4] J. Kapusta and T. Toimela, Phys. Rev. D 39, 3197 (1989). 

[5] E. Braaten and R.D. Pisarski, Nucl. Phys. B337, 569 (1990); B339, 310 (1990). 

[6] J. Frenkel and J.C. Taylor, Nucl. Phys. B334, 199 (1990). 

[7] R.D. Pisarski, Phys. Rev. Lett. 63, 1129 (1989). 

[8] E. Braaten and R.D. Pisarski, Phys. Rev. D 46, 1829 (1992). 

[9] R. Kobes, G. Kunstatter, and K. Mak, Phys. Rev. D 45, 4632 (1992). 
[10] V.V. Lebedev and A.V. Smilga, Physica A181, 187 (1992); T. Altherr, E. Petitgirard, and T. del Ri'o Gaztelurrutia, Phys. 
Rev. D 47, 703 (1993); R. Baier and R. Kobes, ibid. 50, 5944 (1994); K. Takashiba, Int. J. Mod. Phys. All, 2309 (1996); 
J.-P. Blaizot and E. Iancu, Phys. Rev. D 55, 973 (1997); 56, 7877 (1997). 
[11] A. Linde, Phys. Lett. 93B, 327 (1980); D.J. Gross, R.D. Pisarski, and L.G. Yaffe, Rev. Mod. Phys. 53, 43 (1981). 
[12] R.D. Pisarski, Phys. Rev. D 47, 5589 (1993). 



.32 



F. Flechsig, A.K. Rebhan and H. Schulz, Phys. Rev. D 52, 2994 (1995). 

R. Baier, S. Peigne, and D. Schiff, Z. Phys. C 62, 337 (1994); P. Aurenche, T. Becherrawy, and E. Petitgirard, ENSLAPP- 



A-452/93, NSF-ITP-93-155, |hep-ph/940332C . 
F. Flechsig and H. Schulz, Phys. Lett. B349, 504 (1995). 

U. Kraemmer, A.K. Rebhan and H. Schulz, Ann. Phys. (N.Y.) 238, 286 (1995). 
F. Flechsig and A.K. Rebhan, Nucl. Phys. B464, 279 (1996). 

E.V. Shuryak, Zh. Eksp. Teor. Fiz. 74, 408 (1978) [Sov. Phys. JETP 47, 212 (1978)]. 

P. Aurenche, F. Gelis, R. Kobes, and E. Petitgirard, Z. Phys. C75, 315 (1996); A. Niegawa, Phys. Rev. D 56, 1073 (1997); 

P. Aurenche, F. Gelis and H. Zaraket, ibid. 62, 096012 (2000). 

Yu.A. Markov and M.A. Markova, Transp. Theory Stat. Phys. 28, 645 (1999). 

Yu.A. Markov and M.A. Markova, J. Phys. G 26, 1581 (2000). 

J.-P. Blaizot and E. Iancu, Nu cl. Phys. B417 , 608 (1994). 

A. Niegawa, OCU-PHYS-179, frep-ph/0101216 . 

J.C. Taylor and S.M.H. Wong, Nucl. Phys. B346, 115 (1990). 

R. Jackiw and V.P. Nair, Phys. Rev. D 48, 4991 (1993). 

J.-P. Blaizot and E. Iancu, SACLAY-T01/005, CERN-TH/2000-272, [hep-ph/010lTo3 

E. Braaten, R.D. Pisarski, and T.C. Yuan, Phys. Rev. Lett. 64, 2242 (1990); S.M.H. Wong, Z. Phys. C53, 465 (1992). 

F. R. Gantmaher, The Theory of Matrices (Nauka, Moscow, 1988). 
Yu.A. Markov and M.A. Markova, J. Phys. G 27, 1869 (2001). 



